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ABSTRACT 

The Gaussian Hermitian matrix model was recently proposed to have a dual string descrip- 
tion with worldsheets mapping to a sphere target space. The correlators were written as sums 
over holomorphic (Belyi) maps from worldsheets to the two-dimensional sphere, branched over 
three points. We express the matrix model correlators by using the fuzzy sphere construction 
of matrix algebras, which can be interpreted as a string field theory description of the Belyi 
strings. This gives the correlators in terms of trivalent ribbon graphs that represent the cou- 
plings of irreducible representations of su{2), which can be evaluated in terms of 3j and 6j 
symbols. The Gaussian model perturbed by a cubic potential is then recognised as a gener- 
ating function for Ponzano-Regge partition functions for 3-manifolds having the worldsheet as 
boundary, and equipped with boundary data determined by the ribbon graphs. This can be 
viewed as a holographic extension of the Belyi string worldsheets to membrane worldvolumes, 
forming part of a holographic hierarchy linking, via the large N expansion, the zero-dimensional 
QFT of the Matrix model to 2D strings and 3D membranes. 
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A. 4 Summary of the trivalent graph and Pachner moves 



1 Introduction 

The correlators of products of traces in the Gaussian Hermitian matrix model can be expressed 
in terms of triples of permutations that multiply to the identity. This fact, known in the early 
nineties [U |2], was recently revisited in [3]. It was used to propose that the Hermitian matrix 
model has a dual string theory with the 2-sphere S"^ as the target space. The matrix model can 
be viewed as a zero-dimensional quantum field theory, with the Hermitian matrix field living on 
a point. This duality relies on the connection between permutations and branched coverings of 
Riemann surfaces, which leads naturally to an target space. Holomorphic maps branched at 
three points on the sphere are called Belyi maps. They are determined by graphs embedded on 
the covering Riemann surface, also known as ribbon graphs which are related to the double-line 
diagrams of large N expansions. In the mathematics literature on Belyi maps, these graphs are 
also called Grothendieck's dessins d'enfants (HE]. The idea that the worldsheet string theory 
is the standard A-model topological string with sphere target has been developed in [6|, [7j for 
genus zero (planar) worldsheets. Refined counting of these graphs was developed in |8]. 

In this paper, we develop another approach for arguing in favour of S"^ as a target space in 
the dual string theory of the Hermitian matrix model. It is known that the algebra of N x N 
matrices, for any positive integer A^, can be viewed as the algebra generated by matrices Jj, 
representing su{2) in the A^ = 2j ' + 1 dimensional irreducible representation of spin j. This 
is the fuzzy sphere construction [9j, in which the Casimir equation + J2 + = j{j + 1) 
is viewed as a matrix version of the equation defining the sphere through its embedding in 
Euclidean 3-space M^. In this construction, with A^ taken to approach 00, one can recover 
standard field theory actions on the sphere, with an appropriate choice of matrix action. The 
matrix action of interest to us, which is the Gaussian action perturbed by other traces weighted 
by small couplings, can be viewed as a simple topological version of scalar field theory on the 
sphere. Since quantum field theory on the target space of strings is precisely what string field 
theory attempts to construct, we may view the fuzzy sphere construction as providing the string 
field theory for the string theory of Belyi maps. The fuzzy sphere is distinguished among fuzzy 
geometry constructions in that it uses the matrix algebra for any positive integer A^, in contrast 
to fuzzy projective spaces and other fuzzy co-adjoint orbits which use A^ equal to sequences of 
dimensions of representations of higher rank groups (see for example [lO]). 

The fuzzy sphere construction uses fuzzy spherical harmonics which give an su{2) covariant 
basis for the matrix algebra. The matrix model correlators can be expressed in terms of Wigner 
3j and 6j symbols arising in the product of fuzzy spherical harmonics [H]. One of our main 
results is to show that the correlators can in fact be expressed exclusively in terms of sums 
of 6j symbols. This result will be no surprise to readers familiar with spin networks, as our 
method is in fact an adaptation of arguments from the spin network literature. The analogous 
result in the context of spin networks was proved in [12] , and discussed further in [T51 [T3] . 
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It is convenient, for a compact statement of our next result, to restrict attention to cor- 
relators of cubic traces, or equivalently to the Gaussian action perturbed by trX'^. We will 
also restrict attention, in the first instance, to the leading large N limit, where only spherical 
worldsheets contribute. The matrix model computation of these correlators is a sum over pla- 
nar ribbon graphs, which are graphs embedded on the sphere. Each ribbon graph evaluates to 
a power of N. From the fuzzy sphere connection, each ribbon graph can be expressed as a sum 
of 6j symbols with a structure related to the ribbon graph. 

Our second main result is that the sums of 6j symbols, which compute the matrix model 
correlators, can be viewed as partition functions of the Ponzano-Regge model for a 3-manifold 
with a boundary. The 3-manifold is topologically a ball and the boundary is the 2-sphere 
containing the embedded ribbon graph. The Ponzano-Regge state sum model is a model of 
Euclidean gravity in three dimensions which is known to be related to Chern-Simons theory with 
IS0{3) gauge group. It is the g — )■ 1 limit of the Turaev-Viro model, which generates invariants 
of 3-manifolds using sums over representations of the quantum deformation of su{2). The 
computation of the Ponzano-Regge invariant, with the ribbon graph boundary data, chooses 
a cell decomposition of the 3-manifold in terms of tetrahedra which determine the 6j symbols 
being summed. Our construction builds this tetrahedral cell complex (which we call the Belyi 
3-complex) by extending to the 3D bulk a triangulation of the boundary sphere S"^ which is 
well-studied in the context of the Belyi map literature. 

Our third main result is that a membrane extension of a Belyi map can also be found for 
non-planar ribbon graphs. A non-planar ribbon graph can be embedded without intersection 
on a higher genus closed oriented surface, and a triangulation of this surface can be extended to 
a tetrahedral decomposition of a handlebody in three dimensions. We can therefore relate all 
ribbon graphs generated by the Hermitian matrix model to partition functions of the Ponzano- 
Regge model. 

The Gaussian Hermitian matrix model, which can be viewed as a zero dimensional quantum 
field theory, has a two-dimensional dual string theory with 2D worldsheets and 2D target. The 
2D string theory worldsheet can be lifted to a triangulated ball or handlebody in 3D. It is 
therefore appropriate to interpret the 3D space as a holographic lift of the 2D string worldsheet 
to a 3D membrane worldvolume. We thus have a heirarchy of holographies, linking 



OD Matrix model -H- 2D string -H- 3D membrane 



This lifting shares similarities with the construction presented in p!5], since the hologram 
is also a string worldsheet. It is noteworthy that in the context of M-theory there are also 
conjectured hierachies of holographies, which can be viewed by analogy as a precedent for the 
above hierachy. Eleven dimensional M-theory on a 4-torus has a dual which is 5-dimensional 
(0,2) theory [161 [13 [IE]- This in turn has a dual in terms of a large matrix quantum 
mechanics [I9]. Other lower-dimensional formulations of (0,2) theory in 5D and 4D are also 
reviewed in [20]. It would be fascinating to embed the holographic hierarchy of the Gaussian 
Matrix model in a precise manner in M-theory. 

The paper is organised as follows. In the review Section [2] we recall the connection between 
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Hermitian matrix model correlators and permutation triples. We explain how this leads to the 
Belyi map interpretation. Finally we review triangulations associated to Belyi maps, which 
will play an important role subsequently. In particular, for the case where the ribbon graph 
has trivalent vertices, we distinguish two such triangulations, which we call inner and outer in 
anticipation of their roles in the three dimensional picture. The inner Belyi triangulation is 
the dual of the ribbon graph. The outer Belyi triangulation contains the ribbon graph itself, in 
addition to extra vertices and edges added according to specified rules. 

Section |3] reviews the relevant facts about fuzzy spheres and the connection to quantum field 
theory on the 2-sphere. In Section |4] we explain the calculation of correlators in the Gaussian 
matrix model in terms of the fuzzy sphere. This leads to sums involving 3j and 6j symbols. 
We show that the 3js can be summed to give expressions in terms of 6j symbols only. The 6js 
are the basic building blocks of the Ponzano-Regge model. 

In Section [5] we introduce the Ponzano-Regge model and its g-deformed version, the Turaev- 
Viro model, which serves as a regulator. We then explain a prescription for constructing a 
complex, which is a triangulation of the ball , for each planar ribbon graph. The complex is 
built by associating a constituent 3-complex to each vertex of the ribbon graph and gluing these 
constituent complexes together. Since the gluing is determined by the data of the Belyi map, 
we call the complete complex a Belyi 3-complex. The boundary of the complex is the outer 
Belyi triangulation of S'^ associated to the Belyi map, which in particular includes a copy of the 
ribbon graph itself. The interior of the complex contains the inner Belyi triangulation, which is 
the dual of the ribbon graph. We prove, using 6j identities, that the Ponzano-Regge partition 
function of the Belyi 3-complex thus constructed gives the same answer as the contribution 
to the Hermitian matrix integral from the specified ribbon graph. In Section [6] we extend our 
construction of 3-complexes to non-planar ribbon graphs, and prove that the contribution of 
any ribbon graph matches the Ponzano-Regge partition function of the complex constructed 
from the ribbon graph. Section [7] discusses avenues for further research. 

2 Review: The Hermitian matrix model and Belyi maps 

We start by reviewing the Gaussian Hermitian matrix model, following [3] and [21] . The Hermi- 
tian matrix model can be thought of as a quantum field theory in zero space-time dimensions, 
where the observables are correlators of traces of the Hermitian matrix X, invariant under 
X UXW for unitary matrices U. It captures the non-trivial combinatoric structure of 
higher dimensional theories with gauge symmetry, e.g. the gauged Hermitian matrix quantum 
mechanics discussed in [22]. It is also closely related to the combinatorics of the half-BPS sector 
of N = 4 super- Yang Mills theory [231 [21]. We review the description of correlators in terms of 
sums over conjugacy classes of permutation groups, and exhibit some equivalent diagrammatic 
methods of calculating these correlators. We then discuss the string dual of this theory via the 
counting of Belyi maps, and introduce dessins d'enfants as an important tool in visualising this 
duality. 
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2.1 Generating functionals 

We first consider tlie free Gaussian integral over tlie N x N Hermitian matrices 

Me-^*^^', (2.1) 

where — X, and wfiere tfie integral is performed over all the real degrees of freedom of the 
Hermitian matrices, 

N 

DX^Y[dX\ Yl d{ReX'j)d{lmX'.). (2.2) 

k=l l<i<j<N 

As the functional integration is performed over a finite number of variables weighted by an 
exponentially decaying factor, it is wcll-dcfincd and convergent even after insertion of polyno- 
mials in X. It is also invariant under the adjoint action of U(A^), as the action X UXU^ 
preserves the trace of any product of the matrices. 

We can follow the standard procedure for generating functionals of field theories and intro- 
duce a source term J, which is also a Hermitian matrix, 

Z[J] = J Me-^t^^'+t^^''^) (2.3) 

which leads to the propagator 

._ /Me-^trx^XVX^, 

= SlSl (2.4) 
Correlators with more matrix insertions can be calculated using Wick's theorem. For example, 

{X\XIX^%X'-J = S-S^SISI + SISISIS^ + SIS^SISI. (2.5) 
The general form of Wick's theorem here is 

(^1^1 • • • X'tJ = E 61^^51^^ . . . 5- (2.6) 

re [2"] 

where the sum is performed over all the permutations in the permutation group on 2n elements 
that are products of n disjoint 2-cycles. 
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n «2 13 k 



X X X X = X''' X''' X'-'- x'^. 

n n n u 



Jl J2 J3 J4 



Figure 1: X®^ diagram. 
2.2 Combinatoric and diagrammatic methods 

There is a useful method of computing the correlators by thinking of the matrices as linear 
operators with a basis on an A^-dimensional vector space, and representing the operators and 
contractions diagrammatically [3]. 

Consider an A^-dimensional space V with an orthonormal basis |ej). The linear operator 
associated to the matrix X*^- is 

X\e,)=X\\e,). (2.7) 
By extending to the multilinear operator acting on V"®", we have 

(X ® X ® . . . X)|e,, ® e,, ® . . . e,J = X\X% . . . X'l\ej, ® e,, ® . . . e,J, (2.8) 
and by considering the dual vectors, we can write 

(e*i ® e'^ ® . . . e'"|(X ® X ® . . . X)\e,, ® e,, ® . . . e,„) = X\X% . . . X\. (2.9) 

We can introduce a diagrammatic notation for such products of matrices by drawing lines, 
representing the operators, that join labelled points together, representing the vectors and dual 
vectors, as in Figure [TJ We can also denote the contractions as straight lines joining different 
vectors and dual vectors. 



i k i k 




j I 3 I 

Figure 2: Diagram of a propagator and contraction. 
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Figure 3: The general X'^^ correlator, with the indices suppressed. 

The natural gauge invariant operators in this theory are products of traces, as these are 
invariant under the adjoint action of U(A^) on X. Hence, we contract the free labelled indices 
to form correlators of products of traces. For example, 

(trX^) = 6p%6^6l{X\X%X\X'^J, {trXhrX') = S^S%S^S^{X']X%X\X^^J. (2.10) 

The calculation of the correlators can be more clearly performed diagrammatically by adding 
the contractions determining the product of traces to the top of the diagram and identifying 
the upper and lower sides of the diagram. In such diagrams, each loop represents a contraction 
of the form 5* = N. For example, the trX^ operator corresponds to the cyclic contraction of 
four indices, and forms the top half of each of the four diagrams in Figure |4} Wick's theorem 
generates a sum of three diagrams, and as there are three loops in the first and third diagrams, 
and one loop in the second diagram, the correlator (tr(X^)) evaluates to 2N^ + A^. 




Figure 4: The (trX^) correlator. 

The data contained in the different Wick contractions and products of traces can be ex- 
pressed in terms of permutations. We define the action of a permutation a G on V^"' 
as 

a\ej, ® Cj, ® . . . ejj = \ej^^^^ ® e^-^^,, ® . . . ej^^ J. (2.11) 

For example, 

tr(X^) = (e*i ®e*2 ®e'^ ®e**|(X®X®X®X)|ei2 ®ei3 ®ei4 (g)ei,) 
= try®4((X(g)X®X(g)X)(1234)) 

= try®4((1234)(X®X(g)X(g)X)) (2.12) 
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We can therefore write a general multi-trace correlator as a sum over Wick contractions, 

re [2"] 

= (2.13) 

where we have used the abbreviations X = X®^" and tr2„ = trv/®2n, and where is the number 
of disjoint cycles in the permutation 7, e.g. C(i2)(3){4) = 3 for the permutation (12) (3) (4) G S4. 

We define the delta function 6 (a) on a permutation group by setting 5(cr) = 1 if cr is the 
identity permutation and 6{a) = otherwise. Using = C^-i, this allows us to express the 
multi-trace correlator as 

(tr2„(Xa)) =J2 11 5iaT^)N^^. (2.14) 

We note that the correlator is invariant under the action of conjugacy on a, i.e. a — ?■ rar^^ for 
r G S2n- Hence, we could replace a in the delta function with any a' in the conjugacy class of 
a (denoted [a]) and perform the sum over the conjugacy class weighted by its size |[cr]|, 

(tr2n(Xa)) = ^ E E E '5(^W)iV^^. (2.15) 

From the above, we conclude that the observables of this theory have a purely group theoretic 
description, as sums over triples of permutations that multiply to the identity. The combinatoric 
data can be described diagrammatically in several ways. The traditional physics way is to use 
double line diagrams [25] and the closely related ribbon graphs and Grothendieck's dessins 
d'enfants [26]. 

In the double line description, for a contribution to a correlator tr2n(T0"), r G [2"], a G S2n, 
we associate with each disjoint k-cjcle in a a vertex of order k, with each connecting half- 
edge labelled by the numbers in the cycle. We then connect these vertices together with edges 
corresponding to the disjoint 2-cycles in r. The closed loops formed by the double line graphs 
now correspond to the permutation 7 such that ar^ = 1, and hence the evaluation of a ribbon 
graph is N to the power of the number of closed loops in the double line graph. 



tr4 ((13)(24)(1234)) = tr4 ((1432)) 

< — ) 

= N 



Figure 5: The graph associated with the non-planar contraction of tr(X^). 
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The double lines can be shrunk to single lines, and there is no loss of information if we keep 
track of a local cyclic orientation at each vertex. This cyclic orientation can be viewed as being 
derived from an embedding of the graph on a Riemann surface, of the smallest genus that will 
allow the graph to be embedded without intersections. This gives the ribbon graph description 
consisting of vertices and edges, along with cyclic order at the vertices. This leads directly 
to the description as a dessin d'enfant when we subdivide the edges of the ribbon graph by 
introducing a new type of vertex in the middle of each edge. A dessin d'enfant is a bipartite 
graph, i.e a graph with two types of vertices distinguished as black and white with edges only 
linking black to white, that has cyclic order at the vertices. By labelling the edges, we can 
associate a permutation a to the black vertices and a permutation r to the white vertices. In the 
case at hand, the first permutation determines the trace structure and the second permutation 
is a member of the conjugacy class [2"] . Hence the white vertices are always bivalent and the 
graph is called a clean dessin d'enfant. 




Figure 6: The dessin d'enfant equivalent to the above ribbon graph. 

If the connected components of a dessin are drawn without intersection on a collection of 
surfaces of possibly non-zero genus, then the graph partitions the surfaces into distinct faces. 
We can then add to the dessin a third type of vertex, one for each face, that describes the 
permutation 7 such that ^ar — 1. Each face corresponds to a cycle of 7 determined by the 
ordering of the half-edges on its boundary. 

Dessins can also be used to describe Belyi maps, which are holomorphic maps from Riemann 
surfaces onto a sphere branched at three points. Hence, the counting of triples of permutation is 
equivalent to the counting of Belyi maps. In the following section, we review this construction 
and its interpretation as a string theory. 

2.3 Belyi maps 

The Gaussian Hermitian matrix model has a dual string theory in a similar manner to the 
AdS/CFT correspondence. This correspondence is exact, as the combinatoric data of the 
matrix model correlators can be encoded exactly in the branching of holomorphic maps from 
worldsheets to a target space. 
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Consider a surjective holomorphic map / from a Riemann surface S, consisting of a collection 
of connected components of genus Qi, to the complex projective line, or Riemann sphere, P^. 
For a generic point on the sphere, there are 2n preimages on S, where 2n is the degree of the 
map, but for a finite set of points on the sphere there are fewer inverse images. These points on 
the sphere are the branch points of the map, and their preimages on E are the ramification 
points of the map. Now consider a base point on the sphere away from the branch points, 
and draw closed paths from the marked point around each of the branch points. We label the 
preimages of the punctured sphere by the natural numbers {1,2, . . . ,2n}. The holomorphic 
map can be characterised by a permutation in 5'2„ by for each branch point on the target space. 
The permutation is constructed by following the inverse images of the path from the base point 
around the branch point. 

We now specialise to the case where there are three branch points on the sphere. In this 
case, the holomorphic map / is called a Belyi map, and the pair (S, /) is called a Belyi pair. 
If we form a simple loop by combining the loops around each of the three branch points into 
a single loop, then this loop is contractible on the punctured sphere. The preimage of this 
loop must therefore be a collection of 2n disjoint loops, and can be interpreted as the identity 
permutation acting on the set of 2n elements. Therefore, we can characterise the branching of 
a holomorphic map from a Riemann surface to a sphere by a triple of permutations in that 
multiply to the identity permutation, 

0'i(T2cr3 = 1. (2.16) 

This equation still holds after a conjugacy transformation on each of the three elements 
CTj — )■ 70"j7~^, which reflects the arbitrariness of our choice of labelling of the punctured spheres. 
The cycle structure of these permutations is equivalent to the branching profiles about the 
branch points, and hence the Riemann-Hurwitz formula for the covering of a sphere can be 
written as the sum over the Euler characters of the connected components of S, 

^(2 - 2gi) =An- {2n - C^J - {2n - aj - {2n - C^,) (2.17) 

i 

For the case when one permutation is a product of n 2-cycles, the Riemann-Hurwitz formula is 

^(2-2(7,)=C^ + a-n. (2.18) 

i 

We can thus interpret the sums over triples of permutations from the previous section 
as sums over holomorphic maps from a worldsheet to a target space. By using a different 



normalisation of the correlators, we can write (2.15) as 



iV^--"(tr2„(Xa)) = ^ E E E 5(^V7)iV^^+^^-" (2.19) 



2n' \ 'I 2n' 

' o-'e[o-] 7e52n re[2" 
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(2.20) 



/(H,[2"]):S^P1 ' ^1 i 

This notation denotes a sum over maps with branching profiles given by [a] and [2"] , weighted 
by N raised to the power of the Euler characters 2 — 2gi of the connected components of the 
Riemann surface, and where \Aut{f)\ is the order of the group of maps from the Riemann 
surface to itself that satisfy / o = /. 

The structure of the maps from a Riemann surface to the sphere can be visualised by using 
the notion of dessins d'enfants introduced in the previous section. We introduced a dessin as 
representing the data of a triple of permutations that multiply to the identity, but it also has a 
natural interpretation in terms of Belyi maps. Without loss of generality, we can set the branch 
points on the target sphere to be at {0, 1, oo}. We associate the permutation a in the above 
expressions with the branching profile of 0, and the permutation r G [2"] with the branching 
at 1. If the interval [0, 1] on the real line is drawn on the sphere, then the preimage of this 
interval on the Riemann surface produces the dessin associated to the triple of permutations, 
where the black vertices are the preimages of 0, and the white vertices the preimages of 1. 




Figure 7: The dessin d'enfant for the permutations a = (123)(456), r = (16)(25)(34) in S"; 



We can also clarify this mapping and visualise the third permutation by considering the 
preimage of the branch point at infinity. Draw the intervals [l,oo] on the positive real axis 
and [oo, 0] on the negative real axis on the sphere. Denote the preimages of oo by a cross, 
and colour the preimages of the interval [0, 1] in blue, [1, oo] in black, and [oo, 0] in red. The 
real axis partitions the Riemann sphere into two triangles, and the preimages of these triangles 
forms a triangulation of the Riemann surface S. 

The triangles labelled by a plus are preimages of the same triangle on the Riemann sphere, 
and the triangles labelled by a minus are preimages of the other triangle on the Riemann sphere. 
The numbered labelling of the triangle comes from the labels assigned to the preimage of [0, 1]. 
We can read off the permutation 7 = (o"r)~^ from a dessin by writing down the anticlockwise 
cyclic ordering of the triangles labelled with a minus around each preimage of 00. 

Examples of this generated triangulation are given in Figures |9] and 10 for triples of per- 
mutations (cr, r, 7) where cr is a product of 3-cycles. The permutation 7 that satisfies (7x7 = 1 
can be read off by observing the ordering of the minus-labelled triangles around the preimages 
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Figure 8: The image of any Belyi triangulation is a Riemann sphere triangulated by two 
triangles. 



of oo, denoted by crosses. Note that in Figure [T0| the opposite edges on the boundary are 
identified, so there is only one preimage of oo. 




Figure 9: The Belyi triangulation of a sphere generated by the permutations cr = (123) (456), 
r = (16)(25)(34), and 7 = (14)(26)(35). 

We conclude this section by discussing a partitioning of the Belyi triangulation into two 
separate triangulations that will prove useful later. In the main body of this paper, we shall 
only consider connected Ribbon graphs with trivalent vertices, which are equivalent to dessins 
d'enfants specified by one permutation a that is a product of 3-cycles and another permutation 
T that is a product of 2-cycles. 

Consider the white vertices of a Belyi triangulation, associated to a product of 2-cycles r. 
Each white vertex connects to a pair of black edges that connect to cross vertices associated 
to 7, and each white vertex also connects to a pair of blue edges connected to black vertices 
associated to a. Hence, we can remove the white vertices from the diagram by combining their 
connecting pairs of edges of the same colour to generate a new graph with intersecting blue 
and black edges. The new blue edges and black vertices trace out the ribbon graph version of 
the dessin. 

Next, consider only the (new) black edges and their boundary cross vertices. These edges 
and vertices partition the surface of genus g into disjoint contractible faces, one for each black 
vertex in the Belyi triangulation. Since each black vertex is trivalent, and its three connecting 
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Figure 10: The Belyi triangulation of a torus generated by the permutations a = (123) (456), 
r = (14)(25)(36), and 7 = (162435). The opposite edges bounding the rectangle are identified. 




edges intersect each bordering black edge, we conclude that the black edges partition the 
surface into triangles, and hence also form a triangulation of the surface of genus We call 
the triangulation of the surface of genus g generated from the black edges the inner Belyi 
triangulation. This triangulation is the 2D dual of the ribbon graph. 

Finally, consider again the full Belyi triangulation, but with the white vertices and black 
edges removed. This is a graph containing red and blue edges, bounded by black and cross ver- 
tices. Since the removal of white vertices and black edges from the original Belyi triangulation 
is essentially combining pairs of adjacent triangles into new triangles, this process generates a 
new triangulation of the genus g surface. We call this triangulation the outer Belyi triangu- 



lation. We have given two examples of the inner and outer triangulations in Figures 11 and 

m 

In general, any Belyi map generates a full Belyi triangulation of a Riemann surface, but 
the preimages of [1, 00] and of [00, 0] U [0, 1] only form the inner and outer Belyi triangulations 



^Note that if, after removing the white vertices, the graph contains a blue edge connecting to the same black 
vertex at both ends, then the triangulation generated from the black edges will contain faces that resemble cut 
discs. These faces are triangles with two of the edges identified. 
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(a) (b) 
Figure 12: The inner and outer triangulations generated from the Belyi triangulation in Figure 



10 Opposite edges of the boundary rectangle are identified. 



respectively if the branching data at and 1 is specified by a product of 3-cycles and a product 
of 2-cycles. However, this class of Belyi maps will be the ones generated by the correlators of 
the Hermitian matrix model on the fuzzy sphere, and the inner and outer Belyi triangulations 
will eventually lead to an interpretation of this model in three dimensions. 



3 The fuzzy sphere 

In this section we review the relevant necessary facts in the construction of the fuzzy sphere, 
focusing on its description in terms of the fuzzy spherical harmonics, taking definitions and 
results from [SlEIllII]. 

The fuzzy sphere is a family of noncommutative deformations of the algebra of functions on 
the sphere S'^ [9]. The deformation replaces the commuting coordinate functions of the sphere 
with noncommuting operators acting on a vector space. The deformation is performed in a 
manner that preserves the S0{3) symmetry group of the sphere, but loses the notion of a base 
manifold with a continuum of points. 



3.1 Fuzzy spherical harmonics 

The generators of the fuzzy sphere can be defined as operator versions of the coordinate func- 
tions of a 2-sphere embedded in three-dimensional Euclidean space. These generators are the 
N X N matrix generators of the Lie algebra su{2) in an A^-dimension representation, with a 
specific choice of normalisation, defined by the commutation relations 

[xi,Xj] = iXNeijkXk- (3.1) 
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Here, Aat is a noncommutativity parameter defined by 

"1 



4 - 1) 



R[j{j + l)]-^. (3.2) 



In particular, this choice of normahsation of the Lie algebra generators ensures that the 
quadratic Casimir XiXi is equal to R^, and that the coefficient in the commutator tends to 
zero as the dimension of the representation becomes large. We have also introduced the 
half-integer noncommutativity parameter j, related to the dimension of the representation 
by 

N = 2j + 1, jg|. (3.3) 

We will use these two noncommutativity parameters interchangeably throughout. 

The three operators Xi generate the complex matrix algebra of the fuzzy sphere An- Every 
element f E An has a unique expansion 



N-l 

/=E;!/-...^"---^", (3-4) 



1=0 

where the quadratic constraint XiXi = R^ allows the coefficients fi^,,,ii to be taken to be traceless 
and symmetric. This algebra consists of N x N matrices and has dimension A^^, so it is 
equivalent to the algebra of complex N x N matrices. 

An alternative basis of the algebra An is given by the fuzzy spherical harmonics Yim, 
which are deformations of the classical spherical harmonics. The classical spherical harmonics 
are eigenfunctions of the angular momentum operators and J3 and are labelled by their 
eigenvalues, so we expect their fuzzy sphere counterparts to satisfy similar relations. Unlike the 
classical spherical harmonics, there are only finitely many linearly independent eigenfunctions 
of in the fuzzy sphere algebra, and so there are are only finitely many linearly independent 
fuzzy spherical harmonics. 

Denoting the orthonormal basis vectors in a (2j + 1) dimensional representation of su{2) by 
I jyu) , where j is half- integer, we define the fuzzy spherical harmonics by 

yiUj) = V2rTiY,c'ji:jjfi){M, (3.5) 

where Cj^ are real Clebsch-Gordan coefficients, corresponding to the coupling of a pair of 
irreducible representations of su{2) [27\. By using the orthonormality of the \jfi) and some 
standard relations involving sums of Clebsch-Gordan coefficients given in the appendix, it can 
be shown that these operators are eigenfunctions of the angular momentum operators and 
J3 under the adjoint action, that is, they obey 

[Ji,[J,,Yi^]]= adjXadj^{Yim)) = l{l + l)Yi^, 

[J3,yiJ= adj3(y'/„) =mYim. (3.6) 
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We define the inner product on the algebra using the matrix trace, 



if, 9) :=^tr(/t^). 



(3.7) 



Using the Clebsch-Gordan identities (A. 5) and (A. 7) from the Appendix, we can see that 

{yim,yi'm') = Sl,l'Sm,m', (3.8) 

hence the fuzzy spherical harmonics are an orthonormal basis of the algebra with respect to 
this inner product. From the definition (3.5), it can be seen that the Hermitian conjugate of a 
fuzzy spherical harmonic is 



Y,, 



'Im y > ^l-mi 

SO we can state that the trace of a product of fuzzy spherical harmonics is 

tr(^imi^«2m2) = A^(~)™^^«i,«2'^mi -m2- 



(3.9) 



(3.10) 



This relation will be more frequently used than the inner product (3.7) in the following. 
By calculating the product of a pair of fuzzy spherical harmonics, we can write 



V V — \ ^ Ahms Y 
^hmi^ 121712 / J -^/i mi hm-,, ^ hmg 



(3.11) 



'^hmil2'm2 



^ J2 (-)3j+-3^(2/i + 1) (2/2 + 1) (2/3 + 1) ( 

h, 1713 



h h h 
mi m2 —rus 



h h h 
J J j 



(3. 



where the symbols in braces are the Wigner 3j and 6j symbols that describe the coupling of 
irreducible representations of 5^(2) [IT]. More information and background on these symbols 



is given in Appendix A.l 



This construction has shown that the basis of fuzzy spherical harmonics decomposes into a 
sum of irreducible matrix representations of su{2) as 



i©3©5a 

The asymptotic formula for the 6j symbol 



f2N 



lim(-)3^- 



h h h 
J J J 



h h h 




(3.13) 



(3.14) 



shows that the coefficient ^|'Jmi«2m2 fuzzy algebra reproduces the coefficient for the 

classical spherical harmonics algebra in the large N limit. 
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3.2 Quantum field theory on the fuzzy sphere 



Quantum field theories can be constructed on the fuzzy sphere by analogy with those on the 
commutative sphere, see for example [HI |2S1 ES] • We can construct the standard complex scalar 
quantum field theory on the fuzzy sphere by expanding a general member of An as 



Af-l I 



(3.15) 



1=0 m=-l 



As the fuzzy spherical harmonics satisfy (|3.9|), we can set the matrix X to be Hermitian by 

(3.16) 



demanding that the complex conjugate of aim satisfies 



''Im 



ai- 



The partition function is 



(3.17) 



where we integrate over the A^^ real degrees of freedom of the fuzzy sphere with the measure 



N 



1=0 



m=l 



The more commonly considered action for scalar field theory on the fuzzy sphere is 



(3.18) 



(3.19) 



This expression includes a Laplacian, a mass term, and a general potential term, and results 
in the propagator 



jDXe-s 



LU<Jm,m' 



/(/ + l)+/i2- 



(3.20) 



4 The Gaussian Hermitian matrix model as a fuzzy sphere 

The dynamics of scalar field theories on the fuzzy sphere with Laplacians and other terms have 
been considered in several papers. Here, however, we confine ourselves to considering a topo- 
logical theory on the quantum sphere where the Laplacian vanishes. As we could consider the 
terms in a perturbative expansion of the potential to be operator insertions in the correlators, 
we also set the potential to zero, leaving only a mass term. We therefore set /x^ = A^^ to arrive 
at the generating functional 

Z = j DXe-f t^^'. (4.1) 
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This partition function is a Gaussian Hermitian matrix model like the one discussed in Section 
2| but with the A^^ real degrees of freedom rewritten in the su{2) covariant form aim instead 
of the U{N) covariant form X*^-. In the next section, we will prove that the matrix models are 
equivalent. Note that we have chosen a different factor in front of the action for this generating 
function. This will result in different powers of N appearing in the evaluations of the diagrams 
than those of the previous section, but this is just a different choice of normalisation. 



By using (3.10), we see that 



N-1 

tr(X^) = aimO'l'm'tT{YimYi/m' 
l,l'=0 rn=-lrn'=-l' 
N-1 I 

= E «'>'™' (4.2) 



1=0 m= 



and hence calculate the propagator 



j DXe 2 aimtti'm' 



{O'lm.O'l'm') 



e 2 



trx2 



JDX 



We can again retain the U{N) invariance of the original Hermitian matrix model by again 
considering only correlators of products of traces. Using the fuzzy algebra 

^imi^2m2 — ^ ] •^/imi«2m2^3m3 (4-4) 

and the orthonormality of the spherical harmonics, along with Wick's theorem for the variables 

{ai,m, ■ ■ ■ ai,^m2j = Yl n {(^hni^hm,) (4.5) 

t£[2"] disjoint cycles 
(ij) in T 

we can calculate any correlator by writing out explicitly the factors of J^i^^i.^ - ^'^d performing 
the sums over all labels I and m. Alternatively, we can perform the calculations diagrammati- 
cally as in the original Hermitian matrix model. 



4.1 Equivalence of matrix integration and fuzzy sphere path integral 

We introduce a new notation to more clearly exhibit the cyclic structure of the fuzzy spherical 
harmonics, and to simplify the calculations. In Section 3.1, a derivation of the fuzzy sphere 
algebra is given that results in the relation 

^imi^«2m.2 ~ y^ ^Um^h.m-r^h'mzi (4-6) 
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where 



Iimil2m2 



N{- 



\3j+m3 



7(2/1 + 1) (2/2 + 1) (2/3 + 1) 



nil ^2 



h h h 
3 3 3 



(4.7) 



We abbreviate this relation to Y1Y2 = A-^^^Y^, where i = 1, 2, 3 represents a pair of indices 
li,mi, and the repeated upper and lower indices are summed over. The above expression for 
Ai^ can be put into a cyclically symmetric form by defining the lowering and raising operators 



V12 = N{- 
We thus have 



mi 



-m2 5 



12 



N- 



\»Ttl r (T 



'712'? 
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^1 = ^h,h^ 



A123 = = iV2(-)3i7(2/i + l)(2/2 + l)(2/3 + l) 



h h h 
nil ^2 ^3 



h h h 
3 3 3 



(4.8) 



(4.9) 



This now has manifest cyclic symmetry, as it is the trace of three fuzzy spherical harmonics. The 
expressions for the traces in the expansion of correlators can be simplified using this notation. 
For example, the expression for the trace of four spherical harmonics can be written 



tl' (^17711^2 "12 ^37713 



^125^34^! 



(4.10) 



where the repeated 5 represents the sum over the representation labels I5 and representation 
states ms. In addition, a propagator can be written as 



\(^llini ^l2m2 

) = N r] 



1^12 



(4.11) 



and so a general correlator of the model will be a sum over Wick contractions weighted by 
factors of A^, and Aij^. 

We have written the Gaussian Hermitian matrix model variables in terms of fuzzy spherical 



harmonics, but have not derived the form of the Hermitian matrix model measure (2.2) in terms 



of the new variables. To show that the Jacobian for the change of variables is trivial, we will 
show that arbitrary correlators computed in the fuzzy sphere picture, using the standard fuzzy 



sphere measure (3.18), give the same answer as the standard matrix model computation. 



First we recall from the definition (3.5) that the fuzzy spherical harmonics act on the A^- 
dimensional 5^(2) irrep /i = —j, — j + 1, . . . , +j. We abbreviate these vectors to |/i) in the 
following. We next note that a general product of traces can be expressed by a permutation a 
by writing 



tr2„(Xa-) 



(/il/i2 • • •At2n|-^ 
{IJ1H2 ■ ■ ■IJ'2n\X 



Ml v:A'2n 



X\lJ.a{l) 1^(7(2) ■ ■ ■ f^a(2n)) 



(4.12) 
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where the permutation a E S^n acts on the tensor product of vectors in the same way as in the 
Hermitian matrix modeL 

Next, we consider a general contribution to the correlator (tr2„(Xcr)) = (tr2„(crX)). Wick's 
theorem states that 

Te[2"] 

where r = (iiji) . . . {inin) is summed over all products of disjoint 2-cycles, i.e. the ik^jk are all 
distinct integers from 1 to 2n. 

Next, we consider the action of contraction upon a tensor product of spherical harmonics. 



Using the explicit expression (3.5) for the fuzzy spherical harmonics in terms of the basis, 
we write 



AilM2 ^'11^2 



We can simplify this expression using the properties of the Clebsch-Gordan coefficients (A. 6) 



and (A. 7) from the appendix. Hence we calculate 



liTTli 
Illfl2 V\U2 



liirii 

E l^i)(^2| ® |/^2)(/^l|, (4.15) 



A«lM2 

where we have also used the facts that 



'^j^,2 -^>2i-m^~^ ^^-^^i 

and that (j — ^2) is always an integer. We therefore see that the contraction of a pair of indices 
acts like a transposition on the basis vectors. With the complete tensor product of spherical 
harmonics, a contraction of a pair of indices gives 

ri'^Yi ® . . . ® . . . ® y^- ® . . . ® F2n = E ^1 ® • • • ® l^*)^'^^! ® • • • ® I'^i)^'^*! ® • • • ® ^2n- (4.17) 

UiVj 
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We now consider a general contraction r = {hji) ■ ■ ■ {injn), where the numbers ik, jk repre- 
sent distinct integers between 1 and 2n. Using the above relation for all 1 < ik, jk 2n with 
the implicit contraction of the r)'^^ with the Y^, we deduce that 

^hh _ _ _ r]'"HT2n{yi «)...«) Y2nCr) = (/ii . . . fX2nW'^' ■ ■ ■ V'"^"Yi ® . . . (g) Y2n\f^a{l) ■ ■ ■ l^a{2n)) 

tJ-l,---,IJ'2n 

fJ'2n\ {Wi)Wt{1)\ (8) ■ ■ ■ O |j^2n)(j^r(2n)|) |A*a(l) ■ ■ ■ fJ'a{2n)) 

lJ'2n\'T\fJ'a{l) ■ ■ ■ f^a{2n)) 
l^2n\ra\l^l . ..Il2n) 

(4.18) 

This is the same result as for the general correlator in the original Hermitian matrix model. We 
thus see that changing pictures is merely a change of basis, and that changing the variables of 
integration results in a trivial Jacobian. In particular, the ribbon graph and dessins d'enfants 
methods of representing correlator contributions diagrammatically is still valid, where the edges 
now represent pairs of spin labels (/, m). We can therefore conclude that if a and r represent 
a Wick contraction of a correlator, and the associated ribbon graph (or dessin d'enfant) has F 
faces and E edges, then the evaluation of the contribution to the correlator is 

N-^T]''^' . . . r]'^^^tT2EiYi (g)...(g) Y2E0) = N^-^. (4.19) 

A natural normalisation to use for traces of permutations is to weight a correlator described 
by the permutation a e S2E by iV'"'^. For the ribbon graph associated to such a correlator, 
is the number of vertices V. Therefore, using the formula for the Euler character of a closed 
two-dimensional surface, we can write the ribbon graph evaluation 7?. of a connected ribbon 
graph Q in terms of its genus g, 

n[g] := N^-^rf^^^ . . . rf''^HT2E{Yi Y2e(t) = N'^'^^. (4.20) 

Denoting by gi the genus of each connected component of a ribbon graph Q, we can expand a 
general correlator in terms of Ribbon graphs with 

Wick Wick i 

contractions contractions 

In this sum, each Wick contraction determines a ribbon graph Q. The evaluation of the graph 
denoted by TZ[Q], is Y[iN'^~'^^^- In the following, we will speciahse to connected ribbon 
graphs with genus g, which evaluate to iV^~^^. The focus of this paper is on developing a 
three-dimensional interpretation of TZ[Q] using the Ponzano-Regge model. 



fJ.l---fJ.2n 

tJl---fJ2n 
//l.../i2n 

= tr2n(ra). 
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4.2 Trivalent ribbon graphs 



In this section we demonstrate that the contributions to correlators in the fuzzy sphere matrix 
model are naturally understood using trivalent connected ribbon graphs. We have shown 
that there is a correspondence between contributions to correlators in the fuzzy sphere and 
contributions to correlators in the matrix model, as they are both represented by ribbon graphs 
expressing the combinatorial data. For the fuzzy sphere ribbon graphs, each vertex represents 
a trace of fuzzy spherical harmonics, and each edge represents a sum over spin labels weighted 
by a factor of N~^ri^K The traces of general products of fuzzy spherical harmonics are built 
up by the contractions of factors of Aijk with each other using the raising operators rf^, where 
the factor Aijk is a trace of three spherical harmonics. Thus we see that any contribution to a 
correlator can be expressed in terms of products of traces of triples of spherical harmonics. This 
suggests that the ribbon graph of a general correlator has an equivalent ribbon graph where all 
the vertices are trivalent. We show that this indeed the case. 

A general property of ribbon graphs is that a collection of edges around a vertex can be 
'pulled off' in a manner preserving the cyclic ordering at each vertex to generate a new graph 
with an extra edge and vertex, but preserving the number of faces and the genus of the graph. 
Hence the correlator contributions associated to a ribbon graph before and after 'expansion' are 
equivalent. For example, in the expansion shown in Figure 13, a vertex with n outgoing half- 
edges has been expanded into two vertices with 3 and (n — 1) outgoing half-edges respectively, 
connected by a new edge. 




Figure 13: Contraction and expansion of a ribbon graph 



In particular, this expansion can be performed repeatedly until all vertices in the graph 
are trivalent. In such a ribbon graph, all vertices are associated with the factors Aijk-, and all 
edges with the factors Vj'^K We can therefore restrict our attention to correlators of the form 
A^^((trX^)^), which are generated by the expression 



I Z)Xe-f = ^((trX^)^) = iV2-23. (4.22) 



V Wick 

contractions 
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To conclude this section, we consider the example of a ribbon graph in the expansion of 
N{trX^). 





(a) (b) 
Figure 14: A non-planar tr(X^) graph and an expansion of the same graph. 



The ribbon graph in Figure (14(a)) represents the contribution to the correlator 



Ar-2r/^VV^tr(Fiy2l"3W5l6) = A^" = 1. 



(4.23) 



The products of the spherical harmonics in the trace can be written in terms of A^jk symbols 
by using the product rule of the algebra 



(4.24) 



and the choice of which products to take is equivalent to the choice of expansion of the ribbon 
graph. For example, the expansion shown in Figure ( |14(b) ) corresponds to 



N-hT{iY,Y,)iY'Y')iY,Y')) 
N-'A,,^A''^A,'MYaY^Y,) 
N-^A,^''A^^^Aj'^A^ 



(4.25) 



We can develop more understanding of the fuzzy sphere interpretation of the matrix model 
by considering the structure of the sums over these factors Aijk - in particular, their constituent 
Wigner 3j and 6j symbols. 



4.3 Separating the Wigner 3js and 6js 

The sums corresponding to the trivalent graphs are performed over different su{2) represen- 
tations labelled by li and by the states in the representations mj. It is possible to separate 
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the sums out into an mj-dependent part and a /j-dependent part, and to perform the sums 
over the first to arrive at an expression that has no dependence on the states within the 
su{2) representations, but only on the representation labels /j. Although we know that the final 
evaluation of a ribbon graph sum will always be N"^^^^, this decomposition of the sum is still 
a useful approach to take because it results in a link between the Hermitian matrix model and 
theories involving spin networks and the Ponzano-Regge model. 

A general correlator contribution on the fuzzy sphere can be expressed entirely as a sum 
with weights NA123 and N~^ri^'^, where each number i represents a pair of angular momentum 
variables li, 111,1. The contractions between these symbols can be encoded by permutations, or by 
cyclically ordered trivalent graphs, with a factor ^4123 assigned to each vertex and a propagator 
77^^ to each line. This sum could also be written out fully in terms of Wigner 3j and 6j symbols 
with phases and representation dimension weights as follows: recall the definitions 



^123 = A^i(-)'V(2/i + l)(2/2 + l)(2/3 + l) 



mi m2 



^3 

ma 



h h h 
j j j 



(4.26) 



(4.27) 



Now, since we are considering only trivalent graphs with no exterior edges, we know that 
the number of edges E and number of vertices V satisfies 2E — 3V, which means that the total 
factor of N at the front of the expression is ]\iV+^v-2E _ jyV-E^ Hence, we can associate a 
factor of 



N{ 



^1 h h 
mi m2 ma 



h h h 
j j j 



(4.28) 



to each vertex, a factor 



2/1 + 1 



N 



(-)'"^(^Zi,i2<^mi,-r 



(4.29) 



to each half-edge, and perform the sum over all k, rrii. The 2E sets of spin labels (/«, mj) 
correspond to the different half-edges of the trivalent graph, so summing exactly half of these 
labels can immediately reduce the sum to E sets of spin labels, introducing minus signs in the 
3j symbols. This can be represented diagrammatically by assigning orientations to the edges. 
Hence, for a general ribbon graph ^, we arrive at the expression for a correlator contribution. 



n[g] (-)2^-^Ar^-^ n (-)"^'(2^i + 1) n 

lirrii edges 



vertices 



j j j 



. ) , (4.30) 

±mj ±mj ±mfe ' 



where mj appears with a positive sign in the Wigner 3j if the edge is directed towards the 
vertex. 
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We can partition this sum into two parts by considering just the mj-dependent terms, which 
are the 3js and phase factors (— )™, and performing the sums over the labels mj. This expression 
depends purely on the structure of the graph Q and the spin labels k assigned to the edges, 
and is essentially the evaluation of a spin network discussed in [5U]. Hence, we call this part 
of the sum the spin network state sum, or alternatively the 3j sum, associated to a graph 
with /j-labelled edges. The mj-dependent part of the sum for a graph Q is 



rrii edges imv+inoo V 



vertices 



I, 



±mk 



(4.31) 



This sum is invariant under the interchange rrii — )■ —rrii, so the orientations assigned to edges in 
a ribbon graph are arbitrary, and only relevant when constructing and evaluating the expression 



(4.31). We can write the total trivalent ribbon graph evaluation. 



n[g] 



V-E 



J2s[g,{h}]ll{2h + i) n 

li edges vertice 



j j 



h 
j 



(4.32) 



The ribbon graph contribution to a correlator, 7^[^], is here expressed in a factorised form 
containing the spin network state sum. The evaluation of the spin network state sum S[Q] 
cannot usually be performed by inspection, but can be deduced in all cases by employing an 
algorithm of 3j identities in a systematic manner. These 3j identities correspond to operations 
on the trivalent graphs which we call trivalent graph moves, or 3j moves, and will also prove 
to be useful in understanding the 3D interpretation of the graphs in the next section. 

Before presenting the trivalent graph moves, we first discuss two special graphs that are 
relevant to the algorithic evaluation of S[Q]. The simplest trivalent graph with no external 
edges or self-connecting vertices is the 'theta' graph, denoted 9, 



S 




E(-) 



mi,m2,m3 



1. 



mi+m2+m3 



h h h 
nil ^3 ^2 



h h h 

mi 7712 ^3 



h h h 

—mi —1712 —TTls 



h h h 
mi 7712 Tri^ 



(4.33) 



Here we have used only the symmetries and orthogonality properties of the 3js to evaluate 
^[O]. We have also implicitly assumed that the labels (Zi, I2, 13) satisfy a triangle constraint for 
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the 3js to be non-vanishing (see Appendix 



in (4.32) 



A.l), which is enforced by the factor 



^1 h h 
J J J 



The second graph is the tetrahedral network, whose spin network state sum evaluates to a 
6j symbol purely by definition. 



S 




\mi+m2+m3~\-m4~\-m5~\-me I ^1 ^2 ^3 

mi m2 1713 



h 

-mi 



h ^6 

-m5 me 



^3 h h 



h 

-m2 



-me 



k 

-m4 



h h h 
U h h 



(4.34) 



4.4 The trivalent graph moves 

To present the moves on trivalent graphs associated to 3j identities more clearly, we extend 
the definition of 5" to include graphs with external edges. By convention, we do not assign a 
weight of (— )'^' to the external edges, and do not sum over their labels, but will reintroduce 
these required weights and sums when we connect all the external edges together to create a 
complete graph. 

1. The orthogonality relation between two 3js can be expressed as 



S 




4 



m2,m3 



\m2+m3 I h h h 
mi 1712 17l3 



(2/1 + 1)^'^''*'^"'^'"'' 



-m2 



-1714 



-1713 



(4.35) 



By writing the factor (2/i + 1) diagrammatically as a loop, and denoting a delta function 
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(on / and m labels) by a straight line, we can write this expression as 



S 



s 



O 




4 



s 



4 



(4.36) 



We note that Sl^Ti] = 5'[^]5'['H] for a pair of disconnected graphs 0,71. The (—)"'" 
factor that appears in this identity can be interpreted as ensuring that each line in the 
reduced graph has just a single associated factor of (— )™. In addition, by setting li = 1^ 
and introducing the required factor of (— )'"^ we can recover (4.33). We also note that this 



identity cannot reduce down the tetrahedral network from (4.34) to a simpler evaluation. 



since there are no two edges in the graph connected to the same two vertices. 
2. The 3j symbols also satisfy an identity corresponding to the '2-2' move 



E(-)' 



m,3 



^1 h h 

nil ^2 — ""^3 



5^(-)-(2/6 + l)( 



h h h 

7714 nT'S 

h h h 

ms mi mg 



7712 ^4 



-me 



h h h 
h h h 



(4.37) 



This can be expressed diagrammatically as 




3 





o- 



(4.38) 



3. The third 3j identity is associated to the '3-1' move, which reduces three 3j symbols 
to a single 3j and a 6j, 



m4,m^,mfi 



\m4+m5+m6 



/s li Iq 

ms mi —rriQ 



h h h 

rrii ms —m^ 

h h h 

mi 7712 ^3 



Iq I2 U 

tjiq 7712 —rri/i 

h h h 

h h h 



(4.39) 
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(4.40) 



The inverse of this move is called the '1-3' move. 

We note that the 3-1, 2-2 and orthogonality moves are not independent, as the 3-1 move 
could be deduced from the application of the 2-2 move and orthogonality. Alternatively, 
the orthogonality relation could be deduced from the 2-2 and 3-1 moves, provided the 
graph considered has more than two vertices. However, it is useful to include all these 
moves in the set for later applications. 

4. The final necessary 3j move is the 'parity' move, which is the permutation of a pair 
of edges at a vertex. This move is necessary to reduce down non-planar graphs to planar 
graphs, and is not needed for the reduction of planar graphs. 
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(4.42) 



In Section [6} it will also be useful to apply the 2D duals of the first three moves, which are 
called the Alexander moves, or 2D Pachner moves. The 2D dual of a trivalent graph is a 
triangulation of a surface, and as the first three moves do not alter the genus of a graph, the 
Alexander moves relate triangulations of the same surface. It was proved in [21] that any two 
triangulations of a surface of the same genus can be related by a finite series of the Alexander 



moves. These moves are listed in Figure 15 



4.5 Algorithmic evaluation of S[Q] 

The action of these moves on a labelled trivalent graph will generate in the spin network state 
sum evaluation S[Q] a string of factors of {21 + 1), (— )'i+^2+'3^ Q^^d sums over new labels, as 
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Figure 15: The 3-1, 2-2, and orthogonality Alexander moves are dual to the genus-preserving 
trivalent graph moves. 

well as the 6j symbols, which are the evaluations of tetrahedral networks. The tetrahedral 
networks are irreducible under the trivalent graph moves, as any application of the moves on 
them will generate the same Qj or an expression which evaluates to the same 6j. Therefore, 
the evaluation of a general graph by reduction must be a sum of a product of these factors. We 
show algorithmically that all trivalent graphs can be reduced down in this manner. 

A trivalent ribbon graph partitions a surface into vertices, edges, and faces. Each vertex 
is incident to either one, two, or three faces, and each edge is incident to either one or two 
faces. We say that a face is a polygon if it is homeomorphic to a disc when considered with 
its bounding edges and vertices. A necessary and sufficient condition for a face to be a polygon 
is for each edge bounding the face to be incident to two distinct faces. 

The first step in the algorithm is to isolate a polygon of the ribbon graph. Not all ribbon 
graphs possess a polygonal face, but it is always possible to generate such a face from any ribbon 
graph by applying a single parity move. To see this, follow the boundary of a face around until 
a vertex is visited twice, and apply the parity move at this vertex. This will always produce 
a polygon from a non-polygonal face. Since a planar graph always has a polygonal face, a 
planar graph can be reduced without applying the parity move. Also, a non-planar graph will 
eventually reduce down to graph with no polygonal faces, so it is always necessary to apply a 
parity move at least once to evaluate a non-planar graph. Once a polygon has been isolated. 
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we can use a combination of the remaining trivalent moves to remove the polygon from the 
graph and reduce the number of vertices by two. 





Figure 16: The parity move produces a polygon from a face of a non-planar ribbon graph. 



If the polygon is bounded by a single edge, then it is a 'tadpole', and can be removed using 
the 2-2 move and orthogonality relation. Applying the 2-2 move on the edge that connects the 



vertex to a different vertex, as in (4.43) 



S 




Now, using the orthogonality move, and the 6j identity 
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(4.44) 



we can evaluate the sum directly and deduce that 



S 




2/2 + 1 
2/3 + 1 



s 



(4.45) 



If the polygon is bounded by two edges, then it can be removed using the orthogonality 



relation (4.36). If the polygon has three edges, then it can be reduced to a vertex using the 
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3-1 move, as in (4.40). Otherwise, the polygon has four or more edges, and applying the 2-2 



move on adjacent vertices of the face will reduce the number of edges (and vertices) bounding 
the face by one. Performing this move repeatedly will eventually result in a polygon with three 
edges, which can be reduced to a vertex using the 3-1 move. 





Figure 17: The 2-2 move on a polygon reduces the number of edges bounding the polygon by 
one. 



The generation and reduction of a polygon will always reduce the number of vertices of the 
graph by two. Therefore, this procedure will eventually reduce the graph down to a trivial 
loop, which evaluates to a dimension factor (21 + 1). The string of factors and sums that are 
produced in performing these moves gives the final evaluation of the spin network state sum 
S[Q], which is manifestly independent of the labels rrii. In summary, the algorithm is: 

1. Choose a face that is homeomorphic to a disc. If no face is homeomorphic to a disc, then 
apply the parity move to construct such a face. 

2. Apply the 3-1 and 2-2 moves and the orthogonality relation to remove the face. 

3. Repeat these steps until the graph is reduced to a single loop. 



4.6 Examples of 3j and 6j sums 

In this section we have expanded the Gaussian Hermitian matrix model on the fuzzy sphere and 
found that the correlators are described by sums over trivalent graphs. Each trivalent graph 
corresponds to a sum over spin labels rrii) weighted by 3j and 6j symbols, and the sum over 
rriiS weighted by 3js can be evaluated algorithmically. 

We conclude this section by presenting some explicit 3j and 6j sums associated to some 
trivalent graphs. In these cases, we can perform the sum over all the labels using 6j identities 
to confirm the result (4.20). The simplest graph to consider is the theta ribbon graph given in 
Equation (4.33). For this graph, the 3j sum is trivial, and so the total 6j sum is 
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(4.46) 



Orthogonality of the 6js, and the contraint on the range of summation < k < 2j , give the 
expected final answer, 



(4.47) 



Another simple planar graph to consider is the tetrahedral network given in (4.34). For this 
graph, we can state that 
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This 6j sum can be evaluated explicitly by using the Biedenharn-EUiot identity (A. 20), given 



in the appendix, to climate a spin label and a 6j. The remaining sums can be performed by 
using the orthogonality relation. 

Next, we consider the following non-planar graph. 
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We apply the algorithm to find S[Q]. This ribbon graph has no polygonal faces, so we apply a 
parity move on a vertex to deduce that 



s 




(4.51) 



We can now apply the orthogonality relation on this bubble to deduce that 



^^^^ = ^2kTlT^'''' ^ 



(4.52) 



A second parity move on the graph will reduce this trivalent network to a theta network with 
trivial evaluation, hence we deduce that 



(2/4 + 1) 

and hence that the ribbon graph evaluation is 

)i2+«3+i5+;6(2/, + l)(2/2 + l)(2/3 + l)(2/5 + l)(2/6 + l)x 



(4.53) 
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This expression can be evaluated to A^*^ = 1 by using the identity (A. 23) 



5 A three-dimensional interpretation of the Hermitian 
matrix model 

We have expressed every ribbon graph generated by a Hermitian matrix model correlator as a 
sum over representations of su{2) weighted by Wigner 3j and 6j symbols, and as sums over the 
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representation labels ji weighted by 6js and dimension factors. We can find a three-dimensional 
interpretation of these sums by using the Ponzano-Regge model of quantum gravity, which was 
first introduced in [32] and also reviewed in [331 |T3] . In this section, we show that the partition 
functions of the Ponzano-Regge model correspond exactly to these ribbon graph sums when 
certain boundary conditions are applied. 



5.1 The Ponzano-Regge model 



The Ponzano-Regge model is defined by assigning a partition function to any triangulation of 
a 3-manifold, possibly with boundary, with a spin label ji assigned to each edge, and a sum 
performed over all possible values of the spin labels corresponding to internal edges. The sum 
is weighted by the function 
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(5.1) 



hence the partition function 



(5.2) 



interior edges 



is a function of the values of the spin labels on the boundary. 

This partition function takes a very similar form to the sums associated to planar graphs 
in the Hermitian matrix model. By using a judicious choice of labelled cell complex (i.e. a 
labelled triangulation of a manifold), we can reproduce exactly the ribbon graph sums for 
any graph generated by a correlator, which gives us a way of interpreting the zero-dimensional 
combinatoric theory as a three-dimensional topological theory of gravity. We review the features 
of this model necessary to show this correspondence. 




Figure 18: Each 6j symbol in the Ponzano-Regge state sum corresponds to a labelled tetrahe- 
dron. 
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In the Hermitian matrix model, the range of summation of the spins is < k < 2j, where 
li G Z, but in the Ponzano-Regge state sum exhibited above, the spin labels in general range 
over all possible half-integer values up to infinity. However, there are constraints on the ranges 
that these labels can take, which are imposed by the 6j symbols. Recall that the 6j symbol 



Jl J2 33 
Ja J5 J6 



(5.3) 



is zero unless the triples {ji, 72,73}, {ji, is, je}, {j2,j4, je}, 03,74,75} all satisfy the triangle 
inequalities. These triples correspond to the four triangles in the cell complex that border each 
tetrahedroi]|^ This means that if two out of three edges in a triangle are constrained to a finite 
range, then a triangle inequality constrains the third edge to a finite range. Thus, noting that 
the boundary labels are fixed, we see that there is an iterative way of deducing the set of edges 
in a complex whose labels span a finite range. In particular, if all the edges in the complex span 
a finite range, then the state sum must converge, but if any labels are unconstrained then the 
state sum will in general diverge. It can be shown that any triangulation with a vertex in the 
interior must possess spin labels that diverge. The converse statement, that any triangulation 
with no internal vertices must converge, is not true in general, but does hold for all the cases 
that we consider in this paper. 

We wish to construct labelled tetrahedral complexes whose state sums reproduce the ribbon 
graph sums in a systematic manner. While it is possible to do this without introducing internal 
vertices, it is clearer and more systematic to use triangulations with a single internal vertex, and 
to introduce a method of regularising these sums. In the next section we introduce the Turaev- 
Viro partition function, which is a state sum model similar to the Ponzano-Regge model, that 
naturally constrains all the ranges of summation to be finite. We can use this model to recover 
the Ponzano-Regge state sums, and hence the ribbon graph sums, in the 'classical' limit. 



5.2 The Turaev-Viro model as a regulator for the Ponzano-Regge 
model 

The Ponzano-Regge model assigns partition functions to complexes labelled by the irreducible 
representations of su{2). We can deform this model by replacing the representations of su{2) 
in the complex with labelled representations of the quantum deformation of the Lie algebra 
Uq{su{2)), where g is a deformation parameter. The classical algebra is recovered when q is set 
to 1. This deformed algebra has representations analagous to the irreducible representations 
of su{2), labelled by half-integers 7, and containing {2j + 1) states, which can be recoupled 
to generate quantum 37 and quantum 67 symbols |3l]. Unlike su{2), however, the number of 
representations of the quantum algebra is finite whenever g is a root of unity not equal to 1. 

^Note that the tetrahedron associated with a 6j in the Ponzano-Regge model is different from the tetrahedral 
network associated to a 6j in the previous section. The labels associated to edges meeting at a vertex of a 
trivalent graph satisfy a triangle constraint, while the labels associated to a face of a tetrahedron satisfy a 
triangle constraint in the Ponzano-Regge model. These two tetrahedra are dual to each other. 
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Thus, if we demand that g is a root of unity, and replace all the representation-dependent 
expressions in the Ponzano-Regge state sums with their quantum analogues, then the sums over 
representations have finitely many terms and are thus well-defined. This quantity will diverge 
as q tends towards 1 (while still being a root of unity), but there is a natural way of regulating 
this divergence that gives a g-independent limit which coincides with the Ponzano-Regge model 
for the cases where both state sums are convergent. We can thus define the Ponzano-Regge 
model for divergent sums as being the classical limit of the quantum state sum model. 

A more detailed treatment of the quantum state sum model is given in [35]. We present here 
the details of how the relevant quantities, such as summation ranges, representation dimensions, 
and the 6j symbols, deform after being taken to their quantum analogues. 

We take an integer r > 3, and set q 
integer' 



an r root of unity. We define the 'quantum 



\n\ 



g"/2 _ q-n/2 



(5.4) 



that has the property that [n] — )■ n as r — > oo and g — )■ 1, and define the quantum factorials 



H!:= [n][n- I]... [2][l]. 



(5.5) 



We say that a triple of spin labels {ji, 72,73} satisfy the quantum triangle constraints if they 
satisfy the classical triangle constraints with the extra conditions 



J- < (r - 2)/2, 



3i + 32+ k < r - 2. 



(5.6) 



By taking the explicit expression of a 6j symbol in terms of sums and products of factorials 
with triangle constraints given in [361 137] , we can replace the factorials in the definition of a 6j 
symbol with the quantum factorials, and upgrade the triangle constraints to quantum triangle 
constraints, to generate the quantum 6j symbol 



3i 32 33 
3a 3b je 



(5.7) 



This definition coincides with the definition of a quantum 6j given in terms of the recouplings of 
representations of the quantum algebra [38] . The quantum 6j symbol converges to the classical 
6j as g — )■ 1, but crucially is only non-zero for finitely many ji for each value of r. This means 
that if we replace all the 6j symbols in the Ponzano-Regge state sum with quantum 6js, we 
arrive at an always convergent state sum with weight 



n (-)'^i2j. + 1] n 



interior edges 



tetrahedra 



3l 32 J3 
J4 J5 3g 



(5.8) 



Each term in this expression converges to the classical analogue as g — > 1, hence this 
partition function reproduces the original Ponzano-Regge state sum in the g — 1 limit in the 
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cases where the original Ponzano-Regge state sum converges. For the complexes with interior 
vertices, we use the quantum normalisation factor nP' = — 2r/(g^/^ — g"^/^)^ and define the 
Turaev-Viro partition function 

Zg = w^''' J2 W,, (5.9) 

interior edges 

where v is the number of internal vertices in the triangulation. 

By its construction this partition function is finite for any root of unity q, and converges to 
the Ponzano-Regge partition function when the Ponzano-Regge partition function is finite. It 
is more difficult to prove that Zg tends to a finite value for more general complexes, but for the 
classes of manifolds we will be constructing, the g — )■ 1 limit is well-defined, and hence we will 
take this as the definition of the regularised Ponzano-Regge partition function for complexes 
with interior vertices. 

One of the most important properties of the Turaev-Viro and Ponzano-Regge models is 
triangulation independence. Any two triangulations of a 3-manifold that are equal on the 
boundary can be deformed from one to the other by a series of operations on the complex 
called Pachner moves [M]. These moves are mergings and splittings of glued tetrahedra that 
will change the terms that appear in the 6j sums, but due to two identities relating sums 
of products of quantum 6j symbols, these operations will not change the overall value of the 
partition function. 




Figure 19: The 4-1 Pachner Move. 
The identity corresponding to the 4-1 move is 
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Figure 20: The 3-2 Pachner Move. 



and the identity corresponding to the 3-2 move is the Biedenharn-EUiot identity (which also 
holds in the g — )■ 1 limit), 

Vf-)^"f2a + 11 = (5 

Z-^^ ^ ^ 'a in 7a a in i'7 a in iv ?^ ?k iv ?q in 



5.3 Constructing the manifolds associated to planar graphs 

We present a prescription for constructing a 3-complex from a graph Q in such a way that its 



Ponzano-Regge partition function corresponds to the ribbon graph sum (4.32). We denote this 
construction as a mapping A^[^] from graphs to labelled triangulations of manifolds, and the 
partition function Z[A4] of a labelled triangulation as a function from labelled manifolds 
into C. We find that there is a very simple relation between the matrix integral evaluation of 
a graph 7l[Q] and the composition Z o J\A[Q] = Z[Ai[Q]]. Letting V be the number of vertices 
of the graph, we prove that 



N^z[M[g]] = n[g]. 



(5.12) 



In this section we present and discuss the construction of A^[^] for a general planar graph. 
We will see that the construction has the inner Belyi triangulation in its interior and the outer 
Belyi triangulation on its boundary, hence we call A^[^] the Belyi 3-complex of the graph 
Q. We prove that it reproduces the ribbon graph sum TZ[Q] in the next subsection, and discuss 
the construction of manifolds for non-planar ribbon graphs in Section [6] 

We construct the labelled complex A^[^] piecewise from the graph by defining 3-complexes 
associated to each graph vertex, and then gluing the complexes together using the data of the 
graph edges. First, assign an orientation to each edge of the graph. For each vertex, create a 
complex with two tetrahedra glued together on a single face, and then glue on a tetrahedron for 
each half-edge directed towards the vertex. Thus, for a vertex with three incoming half-edges. 
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define 




(5.13) 



and for a vertex with tliree outgoing half-edges, define 




(5.14) 



and similarly for vertices with two or three incoming half-edges. 

In these and all subsequent diagrams of 3-complexes, we use different colours to denote 
different types of spin labels. The red edges will always have the assigned label j, the blue 
edges will have the assigned label 0, and both will always be on the boundary of the complex. 
The black edges in the complex inherit the labels k from the graph edges incident to the graph 
vertex. They will be in the interior of the glued manifold, and the presence of the tetrahedra 
with red j labels will constrain these labels to run over the integers < /j < 2j. Finally, in 
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this section, the green edges in the glued manifold will always be interior edges connected to 
an interior vertex, and their labels, denoted by Roman indices, will run over infinitely many 
values in general. 

We glue together the complexes associated to a pair of graph vertices connected by an edge 
labelled i by identifying the pair of triangles incident to the corresponding labelled edge i. 





Carrying out this gluing procedure for all connections between vertices for the planar graph Q 
will produce a labelled triangulation of the 3-dimensional ball, which will be our definition of 
A^[^]. The black edges will correspond to a triangulation of the sphere dual to the graph Q 
and will lie in the interior of the manifold. The blue edges trace out a copy of the graph Q on 
the boundary of the ball. 

We can see that this combination of red, blue, and black edges in this complex form the Belyi 
triangulations of a sphere discussed in Section |2.3[ Considering the ribbon graph ^ as a dessin 
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d'enfant, we see that the blue and red edges of the 3-complex form the outer Belyi triangulation 
associated to the graph, embedded on the boundary of the complex. The black edges form the 
inner Belyi triangulation associated to the graph, which is embedded in the interior of the 
3-complex. We can therefore view the construction of the labelled complex A^[^] as 'lifting' 
the Belyi triangulation of a sphere into three dimensions. For this reason, we call the labelled 
complex Ai[Q] the Belyi 3-complex associated to a ribbon graph Q. The edges in a complex 
that form the outer Belyi triangulation will always have the same colour-dependent spin labels, 
j for red edges and for blue edges. Hence, in the following text when we refer to an outer 
Belyi triangulation in a complex, we shall also implicitly include the colour-dependent spin 
labellings assigned to the edges. 

5.4 Evaluating the partition function of a Belyi 3-complex 

We evaluate the Ponzano-Regge partition function of a Belyi 3-complex by taking the g — )> 1 
limit of its Turaev-Viro state sum. The term-by-term limit of a sum can be taken if all the 
labels in the summation are constrained to a finite range by the classical triangle constraints. 
Since we shall always take the large r limit, we assume throughout that the noncommutativity 
parameter j <^r. 

The partition functions Zg, and its classical limit Z, are multiplicative. For a pair of disjoint 
labelled complexes A4i and A^2, we have 



This multiplicative rule is modified slightly when the subcomplexes A^i and M.2 are not disjoint, 
but share a boundary. If all the edges on the shared boundary remain on the boundary of the 
glued complex, then the above relation still holds. For a more general gluing, edges that are 
on the boundary of Aii may be in the interior of Aii U A^2; so new weight factors and sums 
need to be introduced. The general gluing procedure for Z is 



where B is the subset of spin labels assigned to lines in dAii fl dAi2 that are in the interior of 
Ail U A^2, and V is the number of vertices that were on the boundary of A^i and M.2 but in 
the interior of the glued manifold. 

This property means that is a sum over the labels associated to the green and 

black edges in the Belyi 3-complex, weighted by 6js. The zero spin label assigned to the blue 
edge and the identity 



Zg[Ml U = Zg[Ml]Zg[M2]. 



(5.17) 



(5.18) 
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(5.19) 



[2j + 1] 
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means that many tetrahedra in the complex have a trivial state sum evaluation, leaving two 
non-trivial tetrahedra associated to each vertex of the ribbon graph. One of these tetrahedra 
generates a 6j with three repeated j labels, but the other tetrahedron, the 'internal' tetrahedron, 
is not as straightforward to interpret. The complex composed of all these interior tetrahedra 
gives a triangulation of a ball bounded by the inner Belyi triangulation, and its associated 
partition function is actually equal to S[Q], the spin network state sum associated to the 
labelled ribbon graph. This is proved in [12] and reviewed in [H], and is only valid for planar 
graphs. 

In this paper, however, we adopt another approach to showing the equivalence of the par- 
tition function of a manifold and its ribbon graph evaluation. We extend the definition of the 
ribbon graph evaluation TZ to include fragments of graphs, and calculate the corresponding 



changes made to Z o J\A and TZ under the trivalent graph moves of Section |4.4| on fragments of 
the graph. We recall the 2-2 move on a pair of vertices. 



6 



(5.20) 



and the 3-1 move 





(5.21) 



These moves will not change the genus of a graph, so we have 



n 



(5.22) 



and 



7^ 




7^ 




(5.23) 



Hence, we wish to show that the corresponding result holds for N'^ Z o A^, which are 



N^ZoM 



(5.24) 



42 



and 



N^ZoM 




NZoM 



(5.25) 



These results can be proved by evaluating the partition function Z on fragments of the graphs 
before and after the moves are applied. This calculation is shown explicitly in Appendix |A.2 



Now we employ the algorithm for reducing down planar graphs that was described in Section 



4.5 For planar graphs with more than two vertices, only the 3-1 and 2-2 moves are required to 
reduce a trivalent graph down to the two-vertex theta graph, so we can perform these moves 
to arrive at 



N'' Z[M[G]] = N'Z[M[Q]]. 



(5.26) 




Figure 21: The gluing of the complex associated to the theta graph 0. 



The 3-complex associated with the two-vertex theta graph is shown in Figure (21), and 
consists of seven tetrahedra, six interior edges, and an interior vertex. Its partition function is 
derived by the limit of the Turaev-Viro partition function, and is 
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(5.27) 

We have omitted here the factors of (— since the 6j triangle constraints force these phase 
factors to be equal to 1. The sums over the edges labelled a, b, and c are unbounded as q 
tends to 1, so we must perform these sums before taking the classical 6j limit. We first use the 
orthogonality relation 
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where the function A is 1 if the spin labels satisfy a triangle constraint, and zero otherwise. 
We next use the relation 



^(_)2a+2fe+2/3[2a + 1][26 + l]A(a,6,/3) = w\-f'[2h + 1], 

ab 

which generates the quantum factor w'^ of the Turaev-Viro model, to deduce that 
N^Z = lim V [2/1 + 1] [2/2 + 1] [2/3 + 1] X 
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(5.30) 



where the previously generated A(Zi, Z2, ^3) has been absorbed into a 6j. All ranges of summation 
are now finite, so we can take the limit term-by-term of this expression, using Equation (5.19) 
to remove the three 'trivial' 6js, to get 
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(5.31) 



This is exactly the same 6j sum given in Section 4^ for the ribbon graph evaluation of 
the theta graph. Hence, we can quote the above to state that the evaluation of this partition 
function is 



N^Z[M[e]] = {2j + lf = N^. 



(5.32) 
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Comparing this partition function with the known evaluation of a ribbon graph of genus zero, 
we can therefore deduce that for any planar graph, we have the relation 



N^Z[M[g]] = 



(5.33) 



As we know that the ribbon graph evaluation of any planar ribbon graph is A^^, we have 
therefore proved that equation (5.12) holds over all planar graphs. 



5.5 Example: The Tetrahedral Graph 

As an example of the equivalence of the Ponzano-Regge state sums and ribbon graph 6j sums, 
we apply the construction to the tetrahedral planar graph. 




This graph has partition function 



N^Z = NUimw-^ y [2/i + l]... [2/6 + !](-; 
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(5.35) 



By using (5.19), the 6js corresponding to the six tetrahedra with a blue line immediately 
evaluate to [A^]~^. The four tetrahedra surrounding the interior vertex correspond to the 6js 
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containing the labels a, b, c, and d. These four tetrahedra can in fact be reduced to a single 
tetrahedron using the 4-1 Pachner move and its associated 6j identity from Section 5^ This 
removes all the dependence on the quantum regularisation factor w'"^ and the spin labels a, b, c, 
and d from the sum, and all the remaining labels are summed over a finite range, so we can 
take the term-by-term classical limit to state that 



N^Z = N-^ ^(2/i + l)(2/2 + l)(2/3 + l)(2/4 + l)(2/5 + l)(2/6 + 1): 
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j j 
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h h h 
J j j 



h h h 
h h h 



. (5.36) 



This is the same 6j sum that appears in Section 4^ for the associated ribbon graph. 

Finally, we note that the Belyi triangulation associated to the ribbon graph is present in this 
complex. The red and blue boundary edges of the complex form the outer Belyi triangulation 
of the sphere, and the black edges form an inner Belyi triangulation of the sphere embedded 
in the interior. If we project the black edges from the interior onto the boundary, we generate 



the full Belyi triangulation shown in Figure 22 




Figure 22: The manifold construction contains the Belyi triangulation. 



6 Higher genus ribbon graphs 



In section 2.3| we defined the inner and outer Belyi triangulations associated with a Belyi map. 
In the previous section, we described a method of generating a 3-complex from a trivalent planar 
ribbon graph, with an outer Belyi triangulation on its boundary and an inner Belyi triangulation 
in its interior. We then showed that it has the same Ponzano-Regge state sum evaluation as its 



ribbon graph sum. In addition, we can see in examples (5.31) and (5.36) that after summing 



out all the labels associated to the green interior lines, the Ponzano-Regge 6j sum takes the 
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same form as the ribbon graph 6j sums (4.46) and (4.48). In general, the Ponzano-Regge 6j 



sum of a Belyi 3-complex will take the same form as the 6j sum associated to a planar ribbon 
graph. This is because the Belyi 3-complex contains a sub-complex, homeomorphic to a ball, 
bounded by the Zj-labelled inner Belyi triangulation, whose Ponzano-Regge partition function 
is equal to S[Q]. 

The aim of this section is to extend the correspondence between Hermitian matrix model 
ribbon graphs and Belyi 3-complexes to non-planar graphs. We show that a labelled 3-complex, 
homeomorphic to a handlebody, with an outer Belyi triangulation on its boundary, can be 
generated from a ribbon graph of arbitrary genus g, and that the Ponzano-Regge partition 
function of this complex and its associated 6j sum satisfy the relation 

Z o M[Q] = N^-^^ = n[g]. (6.1) 

The main obstacle to this generalisation is that the vertex-by- vertex construction of Ai [Q] 
given in Section |5] is no longer valid in the non-planar case. We must therefore present an 
alternative construction of the complex A^[^] for a graph of higher genus, but this construction 
will no longer have the spin network state sum S[Q] appearing explicitly in its partition func- 
tion. However, the non-planar complex A^[^] will possess an outer Belyi triangulation on its 
boundary and an inner Belyi triangulation in their interior, as in the planar case, and we will 
prove that its normalised Ponzano-Regge partition function is equal to its ribbon graph sum. 
We thus consider our alternative construction as the appropriate generalisation of the planar 
Belyi 3-complex construction. 

6.1 Belyi 3-complexes of handlebodies 

The construction of Section |5] associates a complex of tetrahedra to each vertex from a trivalent 
ribbon graph, or equivalently to each black vertex in a dessin d'enfant. The complexes are then 
glued together to construct a triangulation of a ball with an outer Belyi triangulation on its 
boundary, which we call the Belyi 3-complex. We can then use moves on the tetrahedra to 
show that the Ponzano-Regge state sum evaluation of this complex is identical to the 6j sum. 

We run into several problems when we try to apply the same approach to graphs of higher 
genus. Firstly, for even the simplest cases, the complex constructed from the gluing of the com- 



plexes given in Equations (5.13) and (5.14) does not produce the anticipated genus-dependent 



result. For the non-planar graph shown in Figure 23, numerical computation of the Turaev-Viro 
partition function and consideration of 6j asymptotics both show that the partition function 
equals zero in the g — ?■ 1 limit. This cannot be the correct result, as the ribbon graph evaluates 
to N'^~'^3 = 1. In addition, the generated 3-complex is not a topological manifold. The interior 
vertex is a conical singularity, which is not present in the generated triangulations of the ball. 

Due to the triangulation-independence property of the Ponzano-Regge model, once a man- 
ifold and its boundary triangulation are specified, any triangulation of its interior will yield 
the same partition function. The previous construction generates a 3-complex with the correct 
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2 



N/1 



(a) 



(b) 



Figure 23: The simplest non-planar trivalent graph has genus 1, and its dual is a triangulation 
of the torus, where the labelled edges are identified. 

boundary triangulation, but the wrong partition function. Wc can therefore interpret the fail- 
ure of the previous method for higher genus as being caused by the wrong interior manifold 
being generated by the previous construction. 

A natural guess for the correct manifold to consider is a handlebody. A handlebody can 
be formed in three dimensions by embedding a surface of genus g into and considering the 
volume enclosed within the surface. This manifold is compact and has the genus g surface as 
its boundary. We aim to construct from a general genus graph a triangulation of a handlebody 
where the edges of the triangulation on the boundary correspond to the outer Belyi triangula- 
tion. We call any such triangulation of a handlebody with j-labelled red edges and 0-labelled 
blue edges a Belyi 3-complex of the graph, and write M.[Q] to denote a Belyi 3-complex of a 
graph g. 

To create a triangulation of a handlebody with the boundary data given by the ribbon 
graph, we need to use a different method of constructing the complex. For a general higher 
genus graph, we can construct a Belyi 3-complex in three stages. First, find a triangulation of 
the handlebody with the same genus as the graph, and whose triangulation on the boundary 
is dual to the ribbon graph. This boundary triangulation is the inner Belyi triangulation of 
the ribbon graph. In keeping with the conventions established in the previous section, we will 
colour the boundary inner Belyi triangulation edges in black, and the remaining interior edges 
of this complex in green. 

Next, for each triangle on the boundary of the handlebody, glue on a tetrahedron with three 
red edges labelled j. Finally, for each edge of the inner Belyi triangulation, add on a tetrahedron 
with four red j-labcllcd edges and one blue 0-labelled edge, by gluing the tetrahedron on to 
the faces adjacent to the edge of the inner Belyi triangulation. This gluing will move the 
inner Belyi triangulation into the interior of the complex, and hence generate the required 
triangulated handlebody with an outer Belyi triangulation on its boundary. 
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Figure 24: The solid torus can be triangulated with a single tetrahedron by gluing two faces 
after a 27r/3 rotation. 



To exhibit this procedure, we consider the simplest example of a non-planar ribbon graph, 
The graph is of genus one, so we seek a triangulation of the solid torus 



given in Figure 23 



with two triangles on the boundary. In fact, a particularly simple triangulation of the solid 
torus with two boundary triangles is known in the literature on triangulations of 3-manifolds. 
The single tetrahedron triangulation of the solid torus given in |10], displayed here in Figure 
24, is constructed by gluing two triangles together after a twist by 27r/3. In the diagram, the 
triangles labelled by the triples (hjiyh) are identified, and the two distinct triangles labelled 
(^1; h, h) form the boundary of the solid torus. The boundary triangles of this tetrahedron form 
the inner Belyi triangulation of Figure |23) Hence we can add on two tetrahedra with black 



and red edges, and three tetrahedra with black, red and blue edges, to generate the final Belyi 



3- complex shown in Figure 25 



We can check that this complex has the expected Ponzano-Regge evaluation. It is known 
that a genus one ribbon graph evaluates to = 1, so we expect that A^^Z[A^[^]] = 1 for any 
genus one ribbon graph with V vertices. With the appropriate normalisation A^^ added in, the 
Ponzano-Regge partition function of the complex given in Figure 25 is 



(2j + 



1^5^ (2/1 + 1) (2/2 + 1) (2/3 + 1) 



hhh 



h h h 
h h h 



h h h 
j 3 3 



h h h 
3 3 3 



(6.2) 



We can now apply the Biedenharn-Eliot identity (A. 20) on the sum over I3 to get 



(2j + 



- 5^(2/1 + l)(2/2 + l) 

' hh 



h h h 
3 3 3 



3 3 



h 
3 



(6.3) 
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Figure 25: The triangulation of the sohd torus generated by the simplest non-planar graph. 
Three pairs of triangles in the diagram are identified. 



and the orthogonality relation of the 6js (A. 19) reduces this sum to 



The range of summation of li is constrained by A to the integers < /i < 2j, and hence the 
final evaluation of this partition function is 

N^Z = 1, (6.5) 



as required for consistency with (6.1). 



In this example, the triangulation of the solid torus with the correct inner Belyi triangulation 
was simply stated. In general, however, it is hard to find by inspection a triangulation of a 
handlebody whose boundary corresponds to a given inner Belyi triangulation, and it is necessary 
to employ an algorithmic procedure to construct such a complex. 

For a genus one graph, we can construct such a complex by starting from the one-tetrahedron 



triangulation of the solid torus given in Figure 24 and repeatedly adding new tetrahedra onto 



its boundary. This layering of tetrahedra will alter the triangulation on the boundary of the 
complex, but will not alter the topology of the complex. In fact, the layering will alter the 



boundary triangulation by the Alexander moves introduced in Section |4.4[ The fact that it is 
possible to relate any two triangulations of a (2D) torus by a series of Alexander moves means 
that, given any genus one ribbon graph, there exists a sequence of layerings of tetrahedra onto 
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the one-tetrahedron triangulation of the sohd torus that will produce a solid torus triangulation 
with the inner Belyi triangulation associated to the graph on its boundary. 





(c) 

Figure 26: The change to the boundary triangulation from gluing on a tetrahedron across one 
or two, or three faces is dual to a 1-3, 2-2, or 3-1 move. The boundary triangulation is shown in 
black, its dual graph is superimposed in blue, and the other interior edges are drawn in green. 



In Figure 26, we have shown three ways of layering a tetrahedron on to a complex that do 
not change its topology. The black edges denote the boundary of this complex, and the green 
edges lie in the interior of the complex. In addition, we have superimposed the trivalent graph 
dual to the boundary triangulation in blue. We see from the figure that layering a tetrahedron 
on either one, two, or three faces of the complex will alter the boundary triangulation of the 
complex by a 1-3, 2-2, or 3-1 Alexander move. The Alexander moves, introduced at the end of 
Section 4.4, are dual to the trivalent graph moves. A summary of these two-dimensional moves, 
along with the three-dimensional Pachner moves of Section |5.2[ is presented in a diagram in 
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Appendix A.4 

We can now state a procedure for generating a Belyi 3-complex associated to any graph of 
genus one. Given a genus one ribbon graph, we can employ a series of trivalent graph moves 



that will reduce the graph down to the two- vertex graph in Figure 23(a) By considering the 



triangulations dual to this sequence of graphs, and reversing the ordering of the sequence, 
we generate a sequence of triangulations of the 2D torus, each related to its successor by an 
Alexander move. The initial term in this sequence is the triangulation of the torus given 



in Figure 23(b) and the final term is the inner Belyi triangulation of the original ribbon 
graph. Each of the moves in this sequence can be associated with a layering of a tetrahedron 
on to a complex. By starting from the one-tetrahedron triangulation of the solid torus, and 
sequentially carrying out the layering associated to each Alexander move, we generate a series 
of triangulations of the solid torus. The final complex in the sequence will be the required solid 
torus triangulation with the inner Belyi triangulation on its boundary. Finally, we can add 
on the required red and blue labelled tetrahedra to generate a full Belyi 3-complex with the 
topology of a solid torus. 

We demonstrate this procedure by considering the following example of a genus one graph 
with six vertices. This example will be very useful in the next subsection, as it provides a 
means of constructing triangulations of handlebodies of arbitrary genus, and it will be used in 
the evaluation of the partition function of a Belyi 3-complex of general genus. 





(a) (b) 
Figure 27: A six-vertex genus one graph and its dual triangulation of the torus. 



Consider the graph shown in Figure 27(a) Its dual triangulation, shown in Figure 27(b) 



IS 



related to the two-triangle triangulation of the torus by the sequence of Alexander moves shown 



in Figure 28 Hence, starting from the one-tetrahedron triangulation of the torus, and layering 



on a tetrahedron corresponding to each of these moves, we can construct a triangulation of 



the solid torus whose boundary matches Figure 27(b) The complex constructed from these 



layerings is given in Figure 29(a), and its labelled boundary is shown in Figure 29(b) We can 
create a Belyi 3-complex from this complex by adding on a tetrahedron with three red edges 
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for each of the six faces, and a tetrahedron with a blue edge for each of the nine labelled black 
edges on the boundary. This generated Belyi 3-complex contains twenty tetrahedra in total, 
displayed in Figure [30l and has the partition function 




Figure 28: The sequence of Alexander moves that generates the triangulation of the torus in 



Figure 27(b) from Figure 23(b) 





(a) (b) 
Figure 29: The labelled triangulation of the solid torus and its boundary, which is an inner 



Belyi triangulation, that matches Figure 27(b) The interior edges are drawn in green. 
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(6.6) 
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where we have used the identity (5.19) on the nine tetrahedra containing a blue edge. 

This sum can be calculated using the Biedenharn-Eliot and orthogonality relations on the 
6js. For example, considering the 6js containing the label a, we can use 



I4 Ij CL 
I4 I'J Ij 



I4 Ij d 

h k h 



I4 Ij CL 
k b Ig 



h h h 



h h h 
h h h 



(6.7) 



to remove the sum over a and a 6j. This also reduces the number of 6js containing b by one, so 
we can now apply the Biedenharn-Eliot identity (A. 20) on the sum over the label b, to obtain 
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Applying the same identity successively on the labels Iq, Z4, and Is, we generate 
N'Z = 5^ (2/1 + l)(2/2 + l)(2/3 + l)(2/5 + l)(2/7 + l)(2/9 + l)x 
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h h h 
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(6.9) 



Using the orthogonality relation (A. 19) on /s, then It, we obtain 



N'Z = j^j^, J2 (2^1 + + l)(2/3 + l)(2/9 + l)A(/9, J, 



3} 



h h h 
3 3 3 



h h h 
3 3 3 
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(6.10) 



This is now the same sum as in (4.46), so we can state that the final evaluation of the 
partition function of this complex is 



N^Z[M[g]] = 



1. 



(6.11) 



This is the desired result, as it is consistent with the relation (6.1). 



6.2 Evaluating the partition function of any Belyi 3-complex 

In the previous subsection, we gave a method of constructing a triangulation of a solid torus 
from the Belyi triangulation of a graph. We know that ribbon graphs of genus g have the 
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Figure 30: The twenty-tetrahedron Belyi 3-complex associated to the ribbon graph in Figure 
27 , where three pairs of triangles are identified. Its boundary is an outer Belyi triangulation, 
and the black edges trace out an inner Belyi triangulation. 



evaluation TZ[Q] = N"^'^^, and found that the normalised Ponzano-Regge partition function 
iV^Z corresponded with the ribbon graph evaluation in both the examples calculated in the 
previous section. We made the hypothesis that the partition function of a general Belyi- 
triangulated handlebody of genus g is given by N'^ Z = A^^~^^. In this subsection, we first 
discuss how to construct Belyi 3-complexes for a graph of genus greater than one, and then 
prove that the Ponzano-Regge partition function is equal to the ribbon graph evaluation for all 
graphs. 

It is possible to construct a genus g handlebody complex by gluing g solid torus complexes 
together. A genus g handlebody can be created from a solid torus and a genus {g — 1) han- 
dlebody by identifying the two manifolds at contractible discs on their boundaries. Hence it 
is possible to generate a genus g handlebody by identifying g solid tori at contractible discs 
on their boundaries. A triangle on the boundary of a solid torus complex is contractible if 
it is homeomorphic to a disc, or equivalently if it has three distinct vertices. Therefore if we 
take g copies of a triangulation of the solid torus with at least two contractible triangles on its 
boundary, then these complexes can be glued together to produce a triangulation of a genus g 
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handlebody. The complex given in Figure 29 has two contractable triangles on its boundary, 
and hence can be used to construct a handlebody of arbitrary genus. 

From this initial genus g handlebody complex, we can construct a triangulation of the genus 
g handlebody with any genus g inner Belyi triangulation on its boundary. This complex can be 
constructed similarly to the genus one case by using the layerings shown in Figure 26 Any two 
triangulations of a genus g surface can be related by a series of Alexander moves; in particular, 
there exists a sequence of moves that relate the boundary triangulation of the initial complex to 
any inner Belyi triangulation. By sequentially applying the tetrahedral layerings to the initial 
handlebody complex, we will construct a series of complexes whose boundary triangulations 
are related by a series of Alexander moves, and thus can construct handlebody complexes with 
any required inner Belyi triangulation on its boundary. Finally, from this complex we create 
the required Belyi 3-complex by gluing on the labelled red and blue tetrahedra. 

Having discussed how to construct a Belyi 3-complex associated to any ribbon graph, we 
now prove that the Ponzano-Regge partition function of the Belyi 3-complex is determined by 
its genus. We start by showing that the partition functions of any two Belyi 3-complexes for 
graphs of the same genus are equal. We can do this by employing a similar method as was used 
in the planar case, which is by showing that the partition function of A^[^] does not change 
under the application of the 2-2 and 3-1 trivalent graph moves to Q. We do this by taking two 
general Belyi 3-complexes that are almost identical except for at a few tetrahedra, and whose 
boundaries are outer Belyi triangulations of two ribbon graphs related by a trivalent graph 
move. We then show that their partition functions are equal; that is, we show the equalities 
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As the genus-preserving trivalent graph moves are dual to the Alexander moves, it is possible 
to relate any two trivalent graphs of the same genus by a series of the 2-2, 3-1, and 1-3 moves. 
Therefore, for any two graphs, there exists a sequence of Belyi-triangulated complexes A^[^t] 
that differ from their predecessor by a 3-1, 1-3, or a 2-2 move, and hence all have the same 
partition function. Finally, we note that the Ponzano-Regge partition function of a complex 
depends only on the topology of the manifold and the labelled triangulation on its boundary, 
and hence that any two Belyi 3-complexes of the same ribbon graph have the same partition 
function. We can therefore conclude that the partition functions of any two Belyi 3-complexes 
generated from graphs of the same genus are equal. 

As a consequence of this, we now need only calculate the partition function of a single Belyi 
3-complex of each genus to show that any Belyi 3-complex of genus g has the evaluation A^^~^^. 
We can do this by induction on g. First, we take a handlebody Belyi 3-complex of a ribbon 
graph Q of genus g graph with V vertices, and demand that one of the triangles of the inner 
Belyi triangulation is contractible (or equivalently that one of the vertices of the ribbon graph 
lies at the intersection of three distinct faces of the graph). As an inductive hypothesis, we 
assume that its partition function evaluates to iV^~^^. 

Consider the part of the complex near the contractible triangle, which we label {h,l2,h)- 
The complex must consist of a tetrahedron containing the triangle (/i, I2, 13) and three red edges, 
and also three tetrahedra with blue and red edges containing li, I2, and I3, respectively. These 
tetrahedra are shown in Figure 31(a) We can separate the partition function of the complex 
N'^ Z into the partition function of these four tetrahedra and the partition function of the rest 
of the complex by denoting the partition function of the rest of the complex as Z{li,l2,h), so 
we can write its partition function as 



(2j + 1)^ 



^(2/i + l)(2/2 + l)(2/3 + l) 



h h I3 
j j j 



Zik,k,l3) = N^-''- (6.14) 



Note that Z{li,l2,h) is the partition function of a complex formed from the original complex 
by the removal of the four tetrahedra, and has h, I2, and I3 on its boundary. This reduced 



complex is shown in Figure 31(b) 



We can glue the solid torus shown in Figure 29 on to the black triangle of Figure |31(b)| by 
identifying the two triangles labelled {h,l2,h)- Since both triangles are contractible, the new 
complex is topologically the connected sum of a genus g handlebody and a solid torus, and 
hence is a labelled triangulation of a genus {g + 1) handlebody. At this stage, its boundary is 
neither an outer nor inner Belyi triangulation, as its boundary contains five triangles with black 
edges in addition to triangles with blue and red edges. However, we can make the complex a 
Belyi 3-complex by gluing on a collection of red and blue-labelled tetrahedra on to these black 
edges. As before, to each of the black triangles on the boundary we glue on a tetrahedron with 
three red edges (labelled j) and associated 6j 



3 j 3 



(6.15) 
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(a) 



(b) 



Figure 31: The complex near the contractible triangle, and the complex after the removal of 
the four tetrahedra. 



and to each of the black edges on the boundary we glue on a tetrahedra with a blue edge 
(labelled zero), which has the associated 6j 



J J 
k j j 



(2j + l)' 



(6.16) 



In total, we will glue on five tetrahedra that each possess three red edges, and nine tetrahedra 
that each possess a blue edge. The boundary triangulation of the glued complex will be the 
outer Belyi triangulation of a genus {g + 1) ribbon graph Q' with four more vertices than 
the previous ribbon graph Q. As before, we weight the partition function of a Belyi complex 
associated to a ribbon graph with V = V + 4 vertices by a factor of A^^ . We can therefore see 



that its partition function is similar to the expression given in (6.6) for the triangulation of a 



solid torus, but containing a factor of Z, one fewer 6j, and a different power of the prefactor 
A^. The new partition function is 



N^'Z' := AT^+^Z' 
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Zik,k,h)- (6.17) 



Since the partition function Z{li, I2, h) has no dependence on the labels U, h, ■ . ■ , /g, a, b, the 
calculation proceeds identically to in the previous example. Hence, we can follow the previous 



steps (6.6) - (6.10) to state that 



N'^'Z' 
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(6.18) 
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Finally, using the inductive hypothesis (|6.14|), we can state that 



Ar2 



2-29 



N 



2-2(9+1) 



(6.19) 



This is sufficient to prove the conjecture (6.1 ) by induction. Since we have shown in Section 



6.1| that a Belyi 3-complex of genus one with a contractible triangle in its inner Belyi triangu- 
lation has the partition function evaluation A^*^ = 1, we can therefore deduce that any Belyi 
3-complex of any ribbon graph of genus g with V vertices has the partition function 

N^z[M[g]] = N^-^^ = n[g]. 



(6.20) 



7 Discussion 



We now discuss some further aspects of the link between Hermitian matrix models and state 
sum models that we have uncovered, with a view to future research. We start by considering the 
known realisation of the Ponzano-Regge model in terms of Chern-Simons theory with IS0{3) 
gauge group, in the limit of large k, to express our result as a connection between the Hermitian 
matrix model and Chern-Simons. A related line of inquiry would start by using the quantum 
group generalisation of the su{2) algebra to generalise the 3j coupling coefficients that appear 
in our formulation of the Hermitian matrix model. Finally, we consider the non-Hermitian 
complex matrix model, whose correlator contributions are described by bipartite graphs. The 
same methods as described above can be applied to these graphs to generate topological state 
sums, so this suggests a 3D interpretation of the correlators of the half-BPS operators of A/" = 4 
super Yang-Mills in 4D. 



7.1 Euclidean gravity and IS0{3) Chern-Simons theory 



The known links between Chern-Simons theory and the Turaev-Viro model can be used to 
connect our Belyi 3-complex construction to Chern-Simons theory. In this section we briefly 
review these links between Chern-Simons theories, quantum gravity in three dimensions and 
the Turaev-Viro model. More detailed reviews can be found in |^ H2]. 

Quantum gravity in three dimensions was reformulated as a Chern-Simons quantum field 
theory in the eighties [131 HI]. In the frame formalism, the SU{2) one-forms e'* and cu" give a 
Euclidean metric and spin connection on a manifold Ai, and the Einstein-Hilbert action with 
cosmological constant A takes the form 



S 



EH 



M 



e" A {dua + ^eabcW* A u") + ^eabcc" A A e'^ 
/ 6 



(7.1) 



It was shown that this action is equivalent to a Chern-Simons theory with gauge group G, given 
by the action 



Scs[A] 



k 

47r 



tr 



M 



AAdA+-AAAAA 



(7.2) 
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where the gauge group G is 50(4) ~ SU{2) x SU{2) for positive cosmological constant A, 



or I SO (3) when A is zero. It was also shown that the action (7.1) can be rewritten as the 
difference of two Chern-Simons actions with gauge group SU{2), 

Seh = Scs[A+]-Scs[A,], (7.3) 

where 



Al = u''± v^e", A = (^y^ 



(7.4) 



We can therefore see that 3D gravity with A > corresponds to Chern-Simons theory with 
gauge group 5*0(4) and finite level k, and that 3D gravity with A = corresponds to Chern- 
Simons with gauge group 150(3) and — t- oo. After path integration, the partition function 
Zeh associated with the Einstein-Hilbert action with positive cosmological constant and the 
partition functions of 150(3) and SU{2) Chern-Simons theory satisfy the relation 

Zeh - ^cs - \^cs I • V 
and in the large k limit, this becomes 

Z-° = 4f(^)= lim|Z?f)p. (7.6) 

Later, the Reshetikhin-Turaev topological invariant Z^j. for a gauge group G was introduced 
in [15]. This is an invariant of a closed 3-manifold, and is also dependent on a quantum 
deformation parameter k. It is equivalent to the path integral of Chern-Simons theory, which 
can be expressed as Z'^j. = Z^g. 

In [16], the Turaev-Viro partition function was shown to be equal to the absolute square of 
the Reshetikhin-Turaev invariant. Specializing to G = SU{2), 

This can be interpreted as meaning that the Turaev-Viro partition function defines quantised 



gravity in three dimensions given by the action (7.1 ). The equivalence of the Turaev-Viro model 



with 3D quantum gravity was demonstrated explicitly in [17] by showing that the Hilbert space 
of states of the Turaev-Viro model in the A; — )■ oo limit, which is the Ponzano-Regge model. 



is isomorphic to the Hilbert space of states of the theory given by the action (7.1) with A set 
to zero. The space of states in the Ponzano-Regge model is spanned by the labellings of a 
triangulation of a closed 2D surface S, and the Ponzano-Regge partition function for manifolds 
with two S boundary components extends to an inner product on this space. The phase space 
of the 150(3) Chern-Simons theory is the moduli space of flat 150(3) connections on S. This 
is equal to the cotangent bundle of the moduli space of flat SU{2) connections on S. This 
leads to a Hilbert space of the 150(3) Chern-Simons theory spanned by the SU{2) labellings 
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of trivalent graphs of the surface S. An inner product on this space is given by inserting the 
labelled trivalent graphs in a 3-manifold as Wilson operators in the path integral. 

We have thus arrived at an interpretation of the Ponzano-Regge partition functions of our 
Belyi 3-complexes as the partition functions for IS0{3) Chern Simons theory with SU{2)- 
labelled Wilson lines on the boundary. This provides the link between the Gaussian matrix 
model and IS0{3) Chern-Simons. 

There is also an interpretation of the Gaussian matrix model in terms of SU (2) Chern Si- 
mons theory. The correspondence between the Reshetikhin-Turaev and Turaev-Viro partition 
functions of closed manifolds can be extended to manifolds with boundary. In |?T], a general- 
isation of (7.7) is given for a handlebody H, whose labelled boundary triangulation is dual to 
a labelled trivalent graph T. It is stated that 



Z^y^'>[H,T] = Z'^j^^ [Hy^-HJl (7.8) 

i.e. the Turaev-Viro partition function of the handlebody is equal to the Reshetikhin-Turaev 
invariant of the closed manifold H U —H, which is two copies of the handlebody trivially glued 
on their boundary, and with the trivalent graph embedded in the interior. From the point 
of view of the SU{2) Chern-Simons theory, this labelled trivalent graph is a Wilson operator 
embedded in the manifold H U —H. Hence, we can write 

z'cf^'\H, T] = z'cT\h U -H, T]. (7.9) 

In our 3D interpretation of the Hermitian matrix model, Belyi 3-complexes were introduced 
as triangulations of handlebodies with boundary triangulations generated from ribbon graphs. 
From the point of view of the Chern-Simons theory with SU{2) gauge group, this boundary 
triangulation is a network of Wilson lines, each carrying an su{2) representation label j or 0, and 
embedded in a closed 3-manifold. We therefore also have a link between the Hermitian matrix 
model correlators and SU (2) Chern Simons theory on closed 3-manifolds with embedded Wilson 
lines. Chern-Simons terms are known to arise in the worldvolumes of M-theory membranes 
[18] . Our new connections between the Hermitian matrix model and Chern Simons theories 
will hopefully shed light on the embedding of the holographic hierachy of matrix model, Belyi 
maps, and Ponzano-Regge into M-theory. 

There are other known links between theories in zero, two and three dimensions that may be 
usefully compared with our results. For example, Chern-Simons theory with gauge group U{N) 
is related via a matrix model to 2D Yang-Mills theory with U{N) gauge group in |19]. This 
differs from our correspondence, as we have related ?7(iV)-invariant observables in the matrix 
model to IS0{3) (or SU (2)) Wilson operators in 3D. Chern-Simons theory has also been related 
to a matrix model in [50]. Another link between two and three dimensions was observed in 
connection with the perturbative expansion of the Hermitian matrix model in [8] . The refined 
counting of ribbon graphs has a three-dimensional interpretation as a topological field theory 
with permutation groups as gauge groups. Here, our approach using the fuzzy sphere has 
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focused on individual ribbon graphs and has lead to other 3D topological field theories with 
continuous groups {IS0{3) and SU{2)) as gauge groups. Going beyond individual graphs, 
constructing the generating function of the corresponding boundary triangulations using a 
local boundary term in the Chern Simons action is an interesting problem. A more complete 
picture of these diverse 0-2-3 dimensional correspondences would also be highly desirable. 



7.2 A g-deformation of the Hermitian matrix model correlators 

In our correspondence between Belyi triangulations and the Ponzano-Regge model, we used 
the Turaev-Viro model to regulate the divergent Ponzano-Regge state sums. When considering 
triangulations of manifolds with interior vertices, we deformed the 6j symbols representing the 
coupling of representations of su{2) to the quantum 6j symbols representing the coupling of 
representations of Uq{su{2)), where the quantum parameter g is a root of unity. After summing 
out the spin labels that were not bounded from above by the parameter 2j, we took the classical 
limit g — )■ 1 to obtain a finite sum over classical 6js, and hence recovered the Ponzano-Regge 
state sum. The Turaev-Viro model has thus been a mathematical tool to regulate the divergent 
Ponzano-Regge state sum models. 

It is natural to ask if the matrix model correspondence can be extended beyond Ponzano- 
Regge to Turaev-Viro beyond q = 1- For the Belyi complexes of genus g constructed in Sections 
[sjandjoj we have checked that the weighted Turaev-Viro partition function evaluates to [iV]^~^^, 
where the noncommutativity parameter = 2j -|- 1 has been replaced by a quantum integer. 
A matrix model analogue would require an appropriate g-deformation of the Gaussian matrix 
model, which would be defined so as to produce sums over quantum 6js of the type arising 
in Turaev-Viro. One possible prescription would be to use the g-deformed fuzzy spherical 
harmonics described in [51j, but it is not clear that this will lead to an evaluation for ribbon 
graphs of genus g of [A^]^~^^. A successful outcome would lead to a matrix model equivalent of 
finite k Chern-Simons theory with 5*0(4) ~ SU{2) x SU{2) gauge group. 



7.3 Half-BPS operators in the complex matrix model 

In this paper we have developed a correspondence between the Hermitian matrix model and the 
Ponzano-Regge state sum model. We could have instead considered the more general complex 
matrix model with correlator functions 

(tr((TZ®^)tr(rZt®'^)) = j DZ DZ^ e-^*'^^^'tr((TZ®^)tr(rZt®'^), (7.10) 

where the integration is now performed over all A^ x A^ complex matrices, and the permutations 
a,T ^ Sd act on the indices of the matrices. This generalisation doubles the integrated degrees 
of freedom from A^^ real variables to 2A^^ real variables. This model is relevant in A/" = 4 super 
Yang-Mills in 4D, as the matrix model correlators generate the half-EPS operators of a single 
complex scalar field transforming in the adjoint of U(N) [231 152] . 
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The expansion of a general complex matrix in terms of fuzzy spherical harmonics proceeds 
analogously to the Hermitian matrix model, but with a doubling of the degrees of freedom of 
the matrices. The variables and are independent in the complex matrix model, and 
their correlators are 

{aim 0*Vm') = JpSl,l'Sm,m', {aim ai'm') = {olm Oi"'m') = (^-H) 

These relations, together with the separation of traces of Z and Z'^ in the correlator, will cause 
all the graphs generated from the correlators in the complex matrix model to be bipartite. 
As we have not altered the fuzzy spherical harmonics algebra in switching from Hermitian to 
complex matrices, the 3j and 6j sums assigned to each graph remain unchanged. Thus, the 
generalisation to complex matrices results in a more constrained theory where only bipartite 
graphs appear. The 2D-3D correspondence from Belyi triangulations to the Ponzano-Regge 
model still holds for this subfamily of the graphs, so this gives a 3D JS'0(3) Chern-Simons 
interpretation of half-BPS correlators of 4D SYM theory. It would be interesting to see the 
path integal of the Chern-Simons theory arising more directly from the half-BPS sector. 



8 Summary and outlook 

The Gaussian Hermitian matrix model has been a useful toy model for gauge-string duality 
since the nineties [531 iSH EH]. The earliest interpretations recognised it as a model of c = — 2 
matter coupled to Liouville theory, which suggested that the strings exist, in some formal sense, 
in a dual space-time of —2 dimensions. The link between correlators of the Gaussian model and 
Belyi maps was highlighted recently, and used to propose that there is a string interpretation 
with S"^ as target space jS]. As a way to develop the idea of S"^ as a stringy target space, 
we considered the fuzzy sphere construction of matrix algebras, where the matrix becomes a 
quantum field on the fuzzy sphere and the matrix integral becomes a path integral for a scalar 
field on the fuzzy sphere. The computation of correlators in the fuzzy sphere approach leads 
to a labelling of the ribbon graphs of the Hermitian matrix model with su{2) spins, and a sum 
over weights determined by these spins. 

We focused on trivalent graphs which arise in the computation of correlators of (trX'^)^. 
More general correlators can be expressed in terms of these, but the case of cubic graphs allows 
the simplest statement of the key ideas here. The weights of the spin-labelled ribbon graphs 
arising from the multiplication of fuzzy spherical harmonics are given in terms of 3j and 6j 
symbols. We showed that the sums over the 3js can be performed to leave us with sums 
weighted by 6j symbols. We drew on ideas from spin networks, which are closely related to the 
spin-labelled ribbon graphs. 

To describe Belyi maps, it is natural to consider some triangulations of the surface sup- 
porting the ribbon graph that are generated from the graph. In addition to the standard Belyi 
triangulation, we distinguished two other triangulations generated from trivalent ribbon graphs, 
which we called the inner and outer Belyi triangulation. 
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In Section [5} we considered the Belyi triangulations of spheres associated to planar ribbon 
graphs, and extended these to generate triangulations of the three-dimensional ball, with a view 
to describing a precise connection with the Ponzano-Regge model of 3D gravity. The Ponzano- 
Regge model is a state sum model that assigns to a spin-labelled 3-complex a partition function, 
which is a sum over su{2) spins weighted by 6j symbols [32]. The triangulations of the ball are 
three-dimensional complexes of tetrahedra, which we called Belyi 3-complexes, whose edges are 
labelled by representations of su{2). Our construction guarantees that the boundary of the ball 
is equipped with the outer Belyi triangulation of the sphere, while the inner Belyi triangulation 
also appears in the interior of the Belyi 3-complex. Computing the partition function of this 
Belyi 3-complex using identities and sums involving Qjs yields precisely the evaluation of the 
ribbon graph in the Hermitian matrix model. In addition, the Belyi 3-complex possesses a 
sub-complex, also a triangulation of the ball, whose boundary is the inner Belyi triangulation. 
The partition function of this sub-complex corresponds to the evaluation of the 3j sum as a 
product of 6js. Hence the same sum over Qjs is performed in both the Ponzano-Regge and 
Hermitian matrix models. 

The correspondence between the Ponzano-Regge and Hermitian matrix models continues 
to hold for higher genus ribbon graphs, where the triangulation of a higher genus surface is 
extended to a triangulation of a handlebody. The details of this construction are somewhat 
different, as there is no longer a vertex- by- vertex construction of the complex associated to the 
ribbon graph, as there is in the planar case. In addition, the product of 6js corresponding to the 
evaluation of the 3j sum does not appear explicitly in the partition function of a handlebody. 
However, we exhibited the construction of a Belyi 3-complex associated to any ribbon graph, 
which contains the labelled outer Belyi triangulation on its boundary and the inner Belyi 
triangulation, dual to the ribbon graph, in its interior. Finally, we proved that the Hermitian 
matrix model evaluation of any ribbon graph matches up to the partition function of the 
handlebody Belyi 3-complex generated from the graph. 

Several possible extensions of this work have been discussed in Section [7| Given the existing 
conjectures on hierarchies of holographic dualities on M-theory, an embedding of the holographic 
hierarchy of the Gaussian Matrix model into M-theory would be particularly fascinating. 
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A Appendix 



A.l Properties of the Wigner 3j and 6j symbols 

In this appendix we give a summary of the representation theory of su{2) used in this paper, 
taking definitions, expressions and identities from [SHIESIEZ!- The finite dimensional irreducible 
representations of the Lie algebra su{2) are labelled by a half-integer j e {0, |, 1, |, • • •}, where 
the dimension of each representation is (2j + 1). A basis of each representation is \jm), where 
m is a half-integer in {—j, + — The tensor product of two representations 

of su{2) can be decomposed into a direct sum of irreducible representations of su{2). This 
coupling of representations gives rise to the Clebsch-Gordan coefficients to describe how states 
in the coupled representation can be expressed as sums in the uncoupled representation, and 
vice versa. The Wigner 3j symbols are more symmetric versions of these coefficients, and the 
Wigner 6j symbols are a generalisation to describe the coupling of three representations. 

A. 1.1 Clebsch-Gordan coefficients 

For two representations labelled by integer or half- integer ji, j2, a state in the product repre- 
sentation \JM) can be written 

I JM) = J2 nn^Jjimi) ® \j2m,), (A.l) 

mim2 

where the coefficients Cj^^i ^2^2 Clebsch-Gordan coefficients of this representation. 

From constraints such as unitarity and the orthonormality of both representations, and de- 
manding the Clebsch-Gordan coefficients to be real, the inverse equation is 

Ijimi) ® \j2m,) = J^C-nZ nm.\JM). (A.2) 

JM 

The values of J and M in the coupled representation must satisfy 

M = mi+m2, Jg -J2|,|ji -J2I + I,...,ji+J2}, (A.3) 

so the Clebsch-Gordan coefficients are defined to be zero if these conditions are not upheld. 
Note that the conditions 

J3 e {|ji - 32I |ji - j2| + 1, • • • , ji + J2}, ji + J2 + J3 is an integer (A.4) 

are invariant under any permutation of the labels {ji, j2, Js}- As the symmetry group on three 
objects is the triangle symmetry group 5*3, we call these the triangle constraints on the labels 
{ji) J2, js}- We introduce the function A(ji, j2, ja) that is defined to be equal to 1 if the labels 
{ii! J25 js} satisfy the triangle constraints, and zero otherwise. This means that A(ji, j2? Js) = 1 
if the Clebsch-Gordan coefficients C;?^™"^ , ^ are non- vanishing. 

Il IIL~[ 10 III'}, *— ' 
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The coefficients also satisfy the two orthogonahty conditions 



mim2 

and 



J,M 



An exphcit expression for the Clebsch-Gordan coefficients appears in [37J, with which one can 
derive some relations of the coefficients under switching and reversing of their labels. For 
example, most relations can be derived from the three relations 



^jams _ / _\ji+j2-jif]3-m3 

jimi J2m2 V / ji—mi j2—m2 



\ I 12^2 jimi 

fi32-m2 /A Y\ 



2j2 + 1 



A. 1.2 Wigner 3j symbols 



The Wigner 3j symbols are a more symmetric version of the Clebsch-Gordan coupling coeffi- 
cients. We define 



Jl J2 33 

nil 7712 I ^ ^ ^"■'^ ' "^'""^ 



This symbol must still obey the constraint (A.4) to take non-zero values, but now also satisfies 



TUi + 77i2 + 7713 = 0, and has a simplified symmetry under permutation of labels: 
Jl J2 is \ ^ ( is js jl \ ^ ( is jl j2 

7711 7713 J \ ^3 ^1 / \ ^^3 ^li 7712 

l+h+h I 3l js j2 \ _ /_yi+j2+j3 / j2 jl js \ _ (_\h+j2+h I J3 32 Jl 



^ ' mi m3 7712 ) V "^2 f^i iTT'S / ' ' V "^3 "^2 n^i 



\jl+j2+j3 { Jl J2 J3 

-mi — m2 — m3 



(A.9) 



They now satisfy the orthogonality relations 



mim2 

and 



E + 1) ( ,t ) ( ™. 4 ) = (A.IO) 



S + 1) ( ) ( ^, * ) = A(,,. ,3). (A.U) 

jams 
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A. 1.3 Wigner 6j symbols 



The Wigner 3js are associated with the direct sum of two representations with labels ji and 
j2- We now consider the direct sum of three representations. We could first take the sum of 
the representations labelled ji and j2 to make a representation labelled by ji2, and then couple 
on the js to make a triply coupled representation J. Alternatively, we could start by adding 
the j2 and ja to make a j23 representation, and then couple on the ji. 




Figure 32: Two methods of adding the ji, j2, and js representations 



In both cases, we arrive at a triply coupled representation labelled by J. Consider a state 
in the representation \JM) 12 arising from the coupling of ji2 to j^. This can be written 

\JM),2 = J2 '^nfm,2nmjji2mi2)®\j3ms) 

mi2m3 

= E ^r^i. j.m.CjX'lm, \jim,) ® \j2m2) ® \jsms) (A.12) 

The state corresponding to the coupling of J23 to ji is similarly written 

\JM)2s = n^rn2,C=nZ\n.,\ji^i) ® ^2^2) ® ^3^3). (A.13) 

mim2m3m23 

We can consider the difference between the two methods of coupling by considering their inner 
product, 

(TM\ TM\ = \^ f<JM ryji2mi2 f~iJM ^^23^23 /A 1 

\o ivi \i20 ivi 1 2d, / ^ ^ji2mx2 hm^^jxmi j2rn2^jirnx j2-im2z^j2rn2 jzmz \-^-^^) 



mim2m3 
mi2m23 



The state label M has not been summed over in this inner product, but it is shown in [36] that 
this expression is actually independent of M. Therefore, this inner product is a function purely 



67 



of the six j labels. We can now define the Wigner 6j symbol as a symmetrised version of this 
inner product, 



Jl 32 Jl2 
is J j23 



[(2ji2 + l)(2j23 + l)]- 



Vl2+i23 



( JM 1 12 JM) 



23 



[(2ji2 + l)(2j23 + l)] 



\jl2+j2'A 



"112^23 



■yJM 



jimi j2m2^jl2mi2 ^3^3 ^2^2 i3'Tt3 il*^! ^23^23 



E 



\mi+m2+mi2+m:j+M+m2:j 



mi, m2, ms 
mi2, m23, M 



Jl 32 3l2 

mi 7712 fni2 



X 



jl2 J3 J 

—mi2 rris —M 



jl J 323 

—mi M —m23 

J2 J23 J3 

-m2 m23 -ms 



X 



(A.15) 



The 6j symbol is either real or pure imaginary, depending on whether ji2 + 723 is integer or 
half-integer. They are invariant under any permutation of their columns. 
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(A.16) 



and also under the interchange of the upper and lower arguments in any two of their columns. 
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(A.17) 



As they are composed of Clebsch-Gordan coefficients which obey the triangular constraints, 
the above 6j is only non-zero when the sets of variables 



{jl,j2,j3}, {jl,j5,j6}, {j2,j4,j6}, {j3,j4,j5} 



(A.18) 



all obey the constraint (|A.4 ). These triples are mapped into each other by the action of the 
symmetries (A.16) and (A.17). 



A 6j symbol can be viewed as a weight associated with a planar spin network (graph 
embedded on a plane, with edges marked by spins) or as a tetrahedron. The 24 symmetries of 



the 6j generated by (A.16) and (A.17) correspond to the symmetry group of the tetrahedron S^. 



The spin network can be obtained by a 2D duality on the surface of the tetrahedron. In Figure 



33(a) , the trivalent vertices represent the couplings of the triples of representations, and hence 



the triangle constraints are satisfied by the labels meeting at a vertex. This interpretation of 
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the spin network as a 6j is used extensively in Section |4] in describing the 3j sums from ribbon 
graphs. Alternatively, as in Figure 33(b)[ the couplings of representations can be interpreted 
as the triangles of the tetrahedron. This interpretation of the 6j as a tetrahedron is used in 
the Ponzano-Regge model in Section [5] 




(a) (b) 
Figure 33: The 6j symbol as a planar spin network and as a tetrahedron. 



There are many relations and identities regarding sums of 3js and 6js, many of which are 
listed in [56]. Perhaps the most useful Qj relations in this paper are the orthogonality relation 
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A{h,k,k)Mh,h,k)Sj,j„ (A.19) 



and the Biedenharn-Eliot identity, 
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(A.20) 



The Biedenharn-Eliot identity has an interpretation as the equivalence of the Ponzano-Regge 



state sum model under the 3-2 Pachner Move on tetrahedra given in Figure |20 

Also, these relations are useful in explicitly performing calculations of contributions to 
correlators on the fuzzy sphere. 



(A.21) 



J2i2h + 1) 
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A(/2,j,j)A(/3,j,j: 



(A.22) 
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5^(2/i + l) (2/2 + 1) 
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-)2^(2j + l)A(/3,j,j), 



(A.23) 
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(2/i + i)(2/2 + i: 
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(2j + 1)5/3,0, 



(A.24) 
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(2j + 1) 



(A.25) 



A. 2 The planar 2-2 and 3-1 moves in the Ponzano-Regge model 

In this section we demonstrate that the planar moves on ribbon graphs do not change the value 
of the normalised Ponzano-Regge partition function on their associated constructed manifolds. 
We show this locally by considering the complexes associated with the relevant fragments of 
the graph. As in the main text, we freely interchange = 2j + 1 throughout. For the 2-2 
move, we aim to show that 
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we calculate the partition function Z to be 
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and employing the 2-3 move on the two tetrahedra with green lines, we get 



-^$:(2/,+i)(-f«(2/6+i) 



/i I2 h 
h h 



h I3 h 
d b a 



h h h 
d b c 



h I2 h 
J J J 



(A.28) 



Is I4 
j j J 





(A.30) 
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The first of tliese two factors can be manipulated using tlie 3-2 move on its constituent 6js, 
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as required. 

We next consider the 3-1 move 




(A.33) 



and attempt to deduce the corresponding move in terms of Z o M.. Starting from 





(A.34) 
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and using the 3-2 move on the three tetrahedra connected to the green hne labelled d, 



N^ZoM 




X 



2j + 

h h h 
j j J 



J2 (2/4 + l)(2/5 + l)(2/6 + l)x 



/s I4 /s 
J j j 



h h h 
J J J 



h h h 
h h h 



h h h 
a b c 



We perform the sum over the label Iq using the Biedenharn-Eliot identity, 



+ 1) 



h h h 




h h h 


h h I3 




h h I3 




h h h 


3 3 3 




3 3 3 


h h h 




3 3 3 




3 3 3 



and use the orthogonality of the 6js, 



J](2/4 + l)(2/5 + l) 
1415 



h h h 
3 3 3 



h h h 
3 3 3 



to deduce that 



N^ZoM 




{23 + 1)^ 



^1 h h 
3 3 3 



h h I3 
a b c 



(A.35) 



(A.36) 



(A.37) 



Hence 




(A.38) 



(A.39) 
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A. 3 The non-planar 2-2 and 3-1 moves 

In this section, we show that the partition functions of handlebodies generated from ribbon 
graphs related by a trivalent graph move are equaL 

For the 2-2 move, we take two Belyi 3-complexes that are identical except at several tetra- 
hedra on and near their boundary. They have the partition functions 




-Y^{23 + 1) ^(2Zi + l)(2/2 + l)(2/3 + 1) 



x(2/4 + l) (2/5 + 1) 



^1 h h 
J j J 



h h h 
J J J 



z{hMMMM), (A.40) 




Y'{2j + 1) J](2/i + 1)(2Z2 + l)(2/3 + 1) X 



x(2/4 + l) (2/5 + 1) (2/6 + 1) 



I2 U l& 
j j J 



h h h 
J J j 



h h h 
I4 /g 



(A.41) 



where Z(/i, /2, /a, /4, /5) is the partition function of the remainder of the complex not shown 
in the diagram. Since both complete complexes are Belyi 3-complexes, they contain the red 
and blue edges of an outer Belyi triangulation on their boundary. By noting that the blue 
edges trace out two trivalent graphs, we see that the two Belyi triangulated 3-complexes are 
associated to ribbon graphs related by a 2-2 trivalent graph move. 

We can perform the sum over /g in the second partition function using the Biedenharn-Eliot 
identity 



E(2^6 + 1) 



I2 h Iq 
J j j 



h h h 
J J J 



h h h 
h h h 



h h h 
J j j 



h h h 
J J J 



(A.42) 
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and after plugging this identity into |A.41[ we can deduce that the partition functions of the 
two complexes are equal. 

Similarly, for the 3-1 move, we can consider two Belyi 3-complexes that differ only on and 
near their boundary, with partition functions 




'-f' J2 (2/i + l)(2/2 + l)(2/3 + l)(2/4 + 1): 



X (2/5 + 1) (2/6 + 1) 



h h h 
3 3 3 



h h h 
3 3 3 



h h h 
3 3 3 



h h h 
h h h 



Z {h, h, h), 



(A.43) 



NZ 




(2j+l) ^(2/i+l)(2/2+l)(2/3+l) 



h h h 
3 3 3 



Z {h, I2, h)- 



(A.44) 



We can again apply the identity A. 42 to sum out the label Iq from (A.43), resulting in 



N-^Z 




5^ (2/1 + 1) (2/2 + 1) (2/3 + 1) 



X (2/4 + 1) (2/5 + 1) 



h h h 
3 3 3 



h 
3 



h h h 
3 3 3 



Z (/i, h, /s 



h 
3 



h 
3 



(A.45) 
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Finally, we apply the orthogonality relation (A. 23) on the sum over the labels U and I5 to 
deduce that the partition functions are equal. 



A. 4 Summary of the trivalent graph and Pachner moves 



Trivaleiit Graph Moves Pacliuer Moves 




Trivalent graph moves are operations that relate trivalent ribbon graphs on 2D surfaces. 
Genus-preserving trivalent graph moves do not change the genus of the ribbon graph. 
A complete set of genus-preserving trivalent graph moves are sufficient to relate any two 
trivalent graphs of the same genus by a series of moves. An alternative complete set of the 
genus-preserving moves is also given. 
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Pachner moves are operations that relate triangulations of compact manifolds. Any two 
triangulations of a compact manifold can be related to each other by a series of Pachner moves. 

2D Pachner moves are sometimes called Alexander moves, and can relate any two 
triangulations of a compact surface. A complete set of Alexander moves is sufficient to relate 
any two triangulations of a surface of the same genus (oriented surfaces of the same genus are 
homeomorphic). An alternative complete set is also given. The genus-preserving moves on 
ribbon graphs are dual to the Alexander moves. 

The 3D Pachner moves can relate any two triangulations of a compact 3-manifold with 
boundary, provided that the boundary triangulations are identical. A complete set of 3D 
Pachner moves is given that is sufficient to relate any two triangulations of the same 3-manifold. 
An alternative complete set is also given. (The orthogonality 3D Pachner move reduces a 
pair of tetrahedra, glued together on three faces, to a single triangle.) 

References 



[1 
[2: 

[3: 

[4; 

[6: 

[7] 



[9 



M. Bauer and C. Itzykson, "Triangulations", Discrete Mathematics 156, 1996, pp. 29-81. 

E. Looijenga, "Intersection theory on Deligne-Mumford compactifications" , Seminaire Bour- 
baki 35, 1992-1993, pp. 187-212. 

R. d. M. Koch and S. Ramgoolam, "From Matrix Models and Quantum Fields to Hurwitz 
Space and the absolute Galois Group," |arXiv: 1002. 16341 [hep-th]. 

P. Joubert, "Geometric Actions of the Absolute Galois Group", MSc Thesis, 2006. 

L. Schneps, "The Grothendieck theory of dessins d'enfants", Cambridge University Press 
(1994). 



R. Gopakumar, "What is the Simplest Gauge-String Duality?," arXiv: 1104.2386] [hep-th]. 



R. Gopakumar and R. Pius, "Correlators in the Simplest Gauge-String Duality," 
arXiv:1212.1 236l [hep-th]. 



R. de Mello Koch, S. Ramgoolam and C. Wen, "On the refined counting of graphs on 
surfaces," a rXiv: 1209.0334 [hep-th]. 

J. Madore, "The Fuzzy sphere," Class. Quant. Grav. 9 (1992) 69. 



[10] A. P. Balachandran, B. P. Dolan, J. -H. Lee, X. Martin and D. O'Connor, "Fuzzy com- 
plex projective spaces and their star products," J. Geom. Phys. 43 (2002) pp. 184-204 
arXiv:0107099 [hep-th]. 

[11] C. -S. Chu, J. Madore and H. Steinacker, "Scaling limits of the fuzzy sphere at one loop," 
JHEP 0108, 038 (2001) arXiv:0106205v2 [hep-th]. 



78 



[12] J. P. Moussouris, "Quantum models of space-time based on recoupling theory" , D. Phil. 
Thesis (1983). 

[13] R. J. Dowdall, "Spin foam models for 3D quantum geometry", PhD thesis, University of 
Nottingham (2011). 

[14] H. -C. Ruiz, "Toroidal Spin Networks: Towards a Generalization of the Decomposition 
Theorem" , arXiv:1210.0865y l [math-ph] 

[15] A. Dabholkar and S. Murthy, "Fundamental Superstrings as Holograms," JHEP 0802 
(2008) 034 arXiv:0707.3818 [hep-th]. 

[16] T. Banks, W. Fischler, S. H. Shenker and L. Susskind, "M theory as a matrix model: A 
Conjecture," Phys. Rev. D 55 (1997) 5112 arXiv:9610043 [hep-th]. 

[17] N. Seiberg, "Why is the matrix model correct?," Phys. Rev. Lett. 79 (1997) 3577 
arXiv:9710009 [hep-th]. 

[18] A. Sen, "DO-branes on and matrix theory," Adv. Theor. Math. Phys. 2 (1998) 51 
arXiv:9709220 [hep-th]. 

[19] O. Aharony, M. Berkooz and N. Seiberg, "Light-Cone Description of (2,0) Superconformal 
Theories in Six Dimensions", Adv. Theor. Math. Phys. 2 (1998) 119 arXiv:9712117 [hep-th]. 

[20] N. Lambert, C. Papageorgakis and M. Schmidt-Sommerfeld, "Deconstructing (2, 0) Pro- 
posals," arXiv: 1212.3337 [hep-th]. 

[21] P. H. Ginsparg and G. Moore, "Lectures on 2D gravity and 2D string theory" , TASI 1992, 
arXiv:9304011 [hep-th]. 

[22] N. Itzhaki and J. McGreevy, "The Large N harmonic oscillator as a string theory," Phys. 
Rev. D 71 (2005) 025003 arXiv:0408180 [hep-th]. 

[23] S. Corley, A. Jevicki and S. Ramgoolam, "Exact correlators of giant gravitons from dual 
N=4 SYM theory," Adv. Theor. Math. Phys. 5 (2002) 809 arXiv:0111222 [hep-th]. 

[24] D. Berenstein, "A toy model for the AdS / CFT correspondence," JHEP 0407 (2004) 018 
arXiv:0403110 [hep-th]. 

[25] G. 't Hooft, "A Planar Diagram Theory for Strong Interactions," Nucl. Phys. B 72 (1974) 
461. 

[26] Grothendieck, A. (1984), Esquisse d'un programme 

[27] S. Kryszewski and M. Zachcial, "Alternative representation of N x N density matrix", 
arXiv:0602065vl [quant-ph]. 



79 



[28] S. Vaidya, "Perturbative dynamics on the fuzzy and MP^" Phys. Lett. B 512, 403 
(2001) arXiv:0102212vl [hep-th]. 

[29] S. Iso, Y. Kimura, K. Tanaka and K. Wakatsuki, "Noncommutative gauge theory on fuzzy 
sphere from matrix model," NucL Phys. B 604, 121 (2001) arXiv:0101102v4 [hep-th]. 

[30] S. A. Major, "A Spin network primer," Am. J. Phys. 67 (1999) 972 arXiv:9905020 [gr-qc]. 

[31] J. W. Alexander, "The combinatorial theory of complexes", Ann. of Math. 31 (1930) 

[32] G. Ponzano and T. Regge, "Semiclassical limit of Racah coefficients". Spectroscopic and 
Group Theoretical Methods in Physics, North-Holland, Amsterdam (1968) pp. 1-58 

[33] J. W. Barrett and I. Naish-Guzman, "The Ponzano-Regge model," Class. Quant. Grav. 
26 (2009) 155014 arXiv:0803.3319v 2 [gr-qc]. 

[34] H. Ruegg, "A Simple Derivation Of The Quantum Clebsch Gordan Coefficients For 
Sf/(2)g," J. Math. Phys. 31, 1085 (1990). 

[35] V. G. Turaev and O. Y. Viro, "State Sum Invariants of 3-manifolds and Quantum Qj- 
symbols". Topology 31 (1992) pp. 865-902. 

[36] A. R. Edmonds, "Angular Momentum in Quantum Mechanics", Princeton University Press 
(1960). 

[37] M. E. Rose, "Elementary Theory of Angular Momentum", Wiley (1957). 

[38] A. N. Kirillov and N. Y. Reshetikhin, "Representations of the Algebra Ug{sl{2)), q- 
Orthogonal Polynomials and Invariants of Links" , Advanced Series in Mathematical Physics 
7 (1989). 

[39] U. Pachner, "P.L. homeomorphic manifolds are equivalent by elementary shellings", Eu- 
ropean Journal of Combinatorics 12 (1991). 

[40] W. Jaco, J. H. Rubinstein, "Layered-triangulations of 3-manifolds" , arXiv:0603601 [math- 

gt] 

[41] L. Freidel and K. Krasnov, "2-D conformal field theories and holography," J. Math. Phys. 
45 (2004) 2378 arXiv:0205091 [hep-th]. 

[42] S. Carlip, "Quantum gravity in 2+1 dimensions: The Case of a closed universe," Living 
Rev. Rel. 8 (2005) 1 arXiv:0409039 [gr-qc]. 

[43] E. Witten, "Quantum Field Theory and the Jones Polynomial," Commun. Math. Phys. 
121 (1989) 351. 



80 



[44] E. Witten, "(2+l)-Dimensional Gravity as an Exactly Soluble System," Nucl. Phys. B 
311 (1988) 46. 

[45] N. Reshetikhin, V. G. Turaev, "Invariants of 3-manifolds via link polynomials and quantum 
groups", Invent. Math. 103 pp. 547597 (1991). 

[46] V. G. Turaev, "Quantum Invariants of Knots And 3-Manifolds" , de Gruyter Studies in 
Mathematics 18, Waher de Gruyter (1994). 

[47] H. Ooguri, "Partition functions and topology changing amplitudes in the 3-D lattice gravity 
of Ponzano and Regge," Nucl. Phys. B 382, 276 (1992) | hep-th/9112072] . 

[48] O. Aharony, O. Bergman, D. L. Jafferis and J. Maldacena, "N=6 superconformal Chern- 
Simons-matter theories, M2-branes and their gravity duals," JHEP 0810 (2008) 091 
|arXiv:0806.121 8 [hep-th]. 

[49] R. J. Szabo and M. Tierz, "Chern-Simons matrix models, two-dimensional Yang-Mills 
theory and the Sutherland model," J. Phys. A 43, 265401 (2010) |arXiv:1003.1228| [hep-th]. 

[50] M. Mariiio, "Chern-Simons theory, matrix integrals, and perturbative three manifold in- 
variants," Commun. Math. Phys. 253, 25 (2004) arXiv:0207096 [hep-th]. 

[51] H. Grosse, J. Madore and H. Steinacker, "Field Theory on the g-Deformed Fuzzy Sphere 
I", J. Geom. Phys. 38 (2001) 308 arXiv:0005273 [hep-th] 

[5 2] T. W. Brown , "Complex matrix model duahty," Phys. Rev. D 83 (2011) 085002 
|arXiv: 1009.06741 [hep-th]. 

[53] E. Brezin and V. A. Kazakov, "Exactly Solvable Field Theories Of Closed Strings," Phys. 
Lett. B 236 (1990) 144. 

[54] M. R. Douglas and S. H. Shenker, "Strings in Less Than One-Dimension," Nucl. Phys. B 
335 (1990) 635. 

[55] D. J. Gross and A. A. Migdal, "Nonperturbative Two-Dimensional Quantum Gravity," 
Phys. Rev. Lett. 64 (1990) 127. 

[56] D.A. Varshalovich, A.N. Moskalev, and V.K. Khernsonsky, "Quantum Theory Of Angular 
Momentum: Irreducible Tensors, Spherical Harmonics, Vector Coupling Coefficients, 3nj 
Symbols," World Scientific (1989). 



81 



